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ABSTRACT 

Combined gravitational- wave (GW) and electromagnetic (EM) observations of compact binary 
mergers should enable detailed studies of astrophysical processes in the strong-field gravity regime. 
This decade, ground-based GW interferometers promise to routinely detect compact binary mergers. 
Unfortunately, networks of GW interferometers have poor angular resolution on the sky and their EM 
signatures are predicted to be faint. Therefore, a challenging goal will be to unambiguously pinpoint 
the EM counterparts to GW mergers. We perform the first comprehensive end-to-end simulation that 
focuses on: i) GW sky localization, distance measures and volume errors with two compact binary 
populations and four different GW networks, ii) subsequent EM detect ability by a slew of multi- 
wavelength telescopes and, iii) final identification of the merger counterpart amidst a sea of possible 
astrophysical false-positives. First, we find that double neutron star binary mergers can be detected 
out to a maximum distance of 400 Mpc (or 750 Mpc) by three (or five) detector GW networks respec- 
tively. Neutron Star - Black-Hole binary mergers can be detected a factor of 1.5 further out; their 
median to maximum sky localizations are 50-170 deg^ (or 6-65 deg^) for a three (or five) detector 
GW network. Second, by optimizing depth, cadence and sky area, we quantify relative fractions of 
optical counterparts that are detectable by a suite of different aperture-size telescopes across the globe. 
Third, we present five case studies to illustrate the diversity of scenarios in secure identification of 
the EM counterpart. We discuss the case of a typical binary, neither beamed nor nearby, and the 
challenges associated with identifying an EM counterpart both at low and high Galactic latitudes. For 
the first time, we demonstrate how construction of low-latency GW volumes in conjunction with local 
universe galaxy catalogs can help solve the problem of false positives. We conclude with strategies 
that would best prepare us for successfully identifying the elusive EM counterpart to a GW merger. 
Subject headings: binaries: close — catalogs — gravitational waves — stars: gamma-ray burst — 
stars: neutron — surveys 



1. INTRODUCTION 

Our understanding of astrophysical processes in 
strong-field gravity regimes is limited by the nature of 
the electromagnetic (EM) force and current instrument 
sensitivity. Low rates, short timescales, high energies and 
interactions with their environments characterize tran- 
sient strong- field gravity events in the Universe. Such 
events provide us with a brief window to study the inter- 
play of fundamental physical processes. In this respect, 
gravitational wave (GW) astronomy should allow for the 
study of such events in the Universe, currently inacces- 
sible through EM observations. With GW observations 
alone, we may infer the physical and geometric prop- 
erties of individual sources and determine event rates. 
Furthermore, the combination of GW and EM measure- 
ments will lead to improved understanding of astrophys- 
ical processes in the strong-field gravity regime, and to 
also construct a demogr aphic census of different stron^- 
field gravity events (e.g. Bloom et al."2009j Kulkarni & 
|Kasliwal 2009; Phinney 2009). A challenging goal in this 
respect is to localize and identify a strong-field gravity 
event in the Universe jointly by GW and EM measure- 
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ments. 

Within the next decade, a worldwide network of ad- 
vanced versions of ground-based GW i nterferometers — 
LIGO in the US and possibly India ([B arish fc Weiss] 



[19991 ISig g k the TiGO Scientihc Collaboration 2008 
Sathyaprakash & et al. 2012^ , Virgo in Ital y (Xccadia 
,et aL^2011) , KAGRA in Japan ( |Somiya|2012| — will be- 



come operational within the frequency range of 10 Hz to 
a few kHz. At these frequencies, inspiraling and merging 
compact-object binaries, composed of neutron stars (NS) 
and/or stellar-mass black holes (BH), are expected to be 
amongst the most numerous and strongest GW-emitting 
sources. Mergers of NS-NS and NS-BH binaries, where 
at least one object includes initially neutron-rich mate- 
rial and possibly strong magnetic fields, are expected to 
emit in both GWs and EM waves. 

GW detections of compact binary mergers are antici- 
pated to become routine by the end of this decade. Based 
on the observed Galactic binary pulsar distribution and 
population synthesis results, predicted event rates for 
NS-NS binary mergers range from 0.4 to 400 per year 
[with 40 being the mean value quoted in |Abadie et al. 
pOlO)] detectable by an advanced GW three-detector 
network to distances of several hundred Mpc. Based 
solely on population synthesis results due to an absence 
of observed NS-BH systems, predicted event rates range 
from 0.2 to 200 per year for NS-IOM© BH binary mergers 
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with detectable distances > 1 Gpcj^ 

EM detections of compact binary mergers are still a 
matter of debate. There is growing evidence that short 
hard gamma-ray bursts (SGRBs) represent the small 
fraction of NS-NS or NS -BH mergers beamed towards 



us (see e.g. Nakar| |2QQ7| |Berger 2011). Joint GW and 



EM observations can unequivocally test this hypothe 
sis and illuminate the nature of the central engine and 
collimated outflows. Additionally, extensive theoretical 
modeling is underway to predict the EM signature for 
post-merger ejecta, produced from all mergers, that is 
gravitationally unbound to the final remnant BH. For in- 
stance, a plausible EM counterpart in the optical or near 
infrared, referred to as a kilonova, macronova or mini- 
supernova, may be powered by weak ra dioactive decay 
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two scenar- 
ios exist which depend critically on the timescale of the 
EM emission with respect to the NS binary merger time. 
They are: i) a GW event is first detected and is then fol- 
lowed by a slew of mult i- wavelength EM telescopes, and 
ii) an EM observation is seen prior to, or coincident with, 
a GW measurement. Such 'GW-triggered' and 'EM- 
triggered' searches were implemented with the enhanced 
LIGO and Virgo interferometers, before the y halted for 
their u pgrades to their advanced versions (|Abadie et al.| 
2012a[|LIGQ Scientific Collaboration et al.|201:^a|IAbacilel 



et al.[2012b[|LIGO Scientific Collaboration et a l. 2012b). 
In this paper, we locus on the tirst scenario, in which we 
first observe the inspiral of a NS binary in GWs, and we 
then detect their EM counterparts using multiwavelength 
EM observatories. This scenario leverages the instanta- 
neous all-sky visibility of GW detectors in comparison to 
the narrow field-of-view (FoV) of EM facilities. Due to 
the improvement in the instrument sensitivity and hence 
GW-detectable distance by an approximate factor of ten, 
EM follow-up in the advanced GW interferometric era 
presents a new set of challenges to that faced by the ini- 
tial versions. 

Here, we view the EM follow-up of a GW event in three 
steps: i) localization on the sky with GW measured areas 
and volumes using a network of GW interferometers, ii) 
detectability using different multiwavelength EM facili- 
ties, particularly the optical, and iii) strategies to reduce 
the number of false-positive signatures that might mimic 
an EM counterpart of a NS binary merger within the 
same GW localization volume. 

Several works over the past few years have begun to 
explore GW sky localization. A single GW interferom- 
eter has poor directional sensitivity for transient signals 
because of its broad antenna function. Localizing any 
source on the sky depends primarily on triangulating 
the GW signal's arrival times at detectors using net- 
works of three or more GW interferometers. Studies esti- 
mate sky localization errors for unmodeled and modeled 
GW sources to vary from less than one to a few hun- 

The particular value of 10 Mq is chosen to be representative of 
stellar-mass BHs; however, as discussed in § 2, it is unclear whether 
a NS-IOM0 BH binary merger will produce an EM signature. We 
use NS-5M0 BH in this paper. 



dred square degrees using networks of GW interferome- 
ters dFairhurst 2010; Wen & Chen||2010l |Nissanke et al. 

it 



2QTT| 1klimenk Q et al.| |2011| |^chutz||20li| [Veitch et al 
2012[ [Fairhurst|| 2012p . Previous works assume that th 
UW^source is fixed at high signal-to-noise ratio (SNR) or 
distance, and use analytically-de rived Fisher matrix esti- 
mates (e.g.,[Fair hurst 2010; We n fc Chen|2010l|Klimenko 
eFlan[2Qll|lS^utz,,2011, , Veitch et al. | [20 12( |1^ air hurst 
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owever, a significant fraction of expected signals 
at SNR threshold and degeneracies between pa- 
rameters in the predicted GW strain become important. 
In contrast to earlier work, we compute explicit GW er- 
rors of volumes, distances and sky errors using the full 
predicted GW wavestrain for astrophysically distributed 
NS binary mergers. 

Recent works have also begun to explore EM de- 
tectability. Some papers take on a statistical approach; 
[Singer et al.^ (2012) divide sky localization errors be- 
tw een telescopes and advo cate a coordinated response 



or 



Nuttall & Sutton (2010) assign a probabilistic rank- 
ing statistic lor host galaxies out to 100 Mpc . Some 
focus on a p articular w avelength (e.g., X-rays, Kanner 



et al.||2012|). [Metzger fc Bergerl (|2012|) seek to identity 



the most promising amongst proposed EM counterparts 
by defining cardinal virtues and discount foll ow-up in 



optical and radi o compared to 7-rays. There, Metzger 
Bergerl ( [2QT2| assume single numbers for detectable 



distance^ and sky localization. Here, we take on a dif- 
ferent approach. We simulate an astrophysical NS binary 
population and consider the full range of distances, local- 
izations and GW networks. We quantitatively divide the 
pie of binaries by beaming angle. Galactic latitude and 
distance. We then consider the challenges and optimal 
mult i- wavelength strategies in each slice. 

In this work, we present an end-to-end simulation 
with the following five steps jj First, we construct 
astrophysically-motivated distributions of NS-NS and 
NS-5 M0 BH binary mergers detectable by different GW 
networks using different triggering criteria (§2). Each bi- 
nary will have specific geometric properties: an orienta- 
tion, a sky position and a luminosity distance. Second, 
by simulating GW data streams using analytically mod- 
eled GW strains, we estimate source parameters mea- 
sured by different GW networks. For parameter estima- 
tion, we use MCMC methods developed in Nissanke et al.| 
[2010) [2011 (henceforth, NIO and Nil respectively). 7^ 
well as estimating the sources' sky areas, we compute the 
sky volume errors (§3). We summarize the distributions 
of sky errors, volumes and distances for NS binary merger 

^ For clarity, we define here different distance definitions that 
are used in the hterature. By average detectable distance, we refer 
to the average distance that a single GW interferometer with ide- 
alized Gaussian instrument noise can observe NS binary inspirals 
averaged over all possible sky positions and binary orientations. 
On the other hand, the horizon distance refers to the maximum 
detectable distance that a single GW interferometer with idealized 
Gaussian noise can detect a NS binary event that is located directly 
above the interferometer and is optimally orientated face-on. The 
horizon distance improves on the average detectable val ue bv an 
approximate geometric factor of ~ 2.24 (e.g. see Finn &: Chernoffl 
[1993]). 

^Discussed in detail throughout the paper, our results are neces- 
sarily limited by the assumptions we make; for instance, we assume 
that joint observable GW-EM events are non-spinning NS-NS and 
NS -5 M0 binary systems, and idealized Gaussian GW interfero- 
metric noise and observing conditions at optical telescopes. 
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populations detected by different networks and different 
trigger criteria (§4). Third, armed with locahzations and 
distances of each binary in the simulation, we assess the 
feasibility of detecting an EM counterpart with a wide 
suite of current and planned EM facilities. We pay close 
attention to the trade-off between depth and area given 
finite telescope time (§5). Fourth, using detailed case 
studies, we present the challenges and discuss possible 
strategies to pinpoint the GW event amongst the antic- 
ipated few to many false-positive transients in different 
wavelengths (§6). In conjunction with a galaxy cata- 
log, we discuss how fractional reductions in volume error 
can reduce the overall number of false-positive transients 
within a GW- localized event. Finally, we conclude with 
strategies that maximize the success in identifying EM 
counterparts to GW events (§7). 

2. CONSTRUCTION OF TWO NS BINARY MERGER 
CATALOGS 

We construct two distinct catalogs with either 4 x 10^ 
NS-NS or 3 X 10"^ NS-5 BH binary populations^ For 
every binary in each catalog, we assign physical and ge- 
ometric source properties as described below. In this 
study, physical source parameters are the individual com- 
pact objects' masses mi and m2 (the spin of the NS 
and/or BH is assumed negligible). The geometric source 
parameters comprise the luminosity distance Dl^ the 
sky position in spherical polar coordinates (which points 
from the center of the Earth to the binary) n = (^,^), 
and the binary's inclination angle cos ^ = L • n, where L 
is the unit vector normal to the binary's orbital plane. 
The colatitude and longitude (j) describe n, and are 
related to the declination S and right ascension a, by 
= 7r/2-S and = a-GAST respectively, where GAST 
is Greenwich Apparent Sidereal Time (see NIO for details 
on the binary's and Earth coordinate systems used in this 
work) . Let us consider now how we assign specific source 
parameters to each binary. 

Regarding the physical source parameters, we assume 
that each NS has a physical mass of 1.4 Mq and each 
BH has a physical mass of 5.0 M©, and that the ob- 
jects are non-spinning. In practice, we expect the NS 
binary population in the Universe to have continuous NS 
and/or BH mass distributions. In the case of NS-BH bi- 
naries, instead of the fiducial 10 Mq BH used in standard 
GW literature, we choose BHs with a small enough mass 
that the NS companion does not plunge directly into the 
gravitational potential well of the central BH. Therefore, 
we can expect some tidal disruption of the NS to oc- 
cur and to observe an accompanying EM counterpart. 
Tidal disruption occurs if the tidal disruption radius is 
greater than the BH's innermost circular orbit (ICO); the 
tidal disruption radius being a function of the BH's spin, 
the NS's equation of state and the binary's mass ratio 
(see discussions in e.g. [Taniguchi et al . 2007; Shibata & 
Taniguchi|2QQ8l|Shibata et al |:^()()()1 |!^utoku et al.,2011 



6GM/c^ for a non-spinning BH, where G is the Gravi- 
tational Constant, c is the speed of li^ ht and M is the 



BH's mass (jShapiro & Teukolsky 



Foucart et al.||2011| |Foucart||2012p . The ICO describes 
the last stable circular orbit of the binary system prior 
to the merger and can be approximated to a test par- 
ticle's innermost stable circular orbit (ISCO) radius of 

^ We choose a sufficiently large number of binary systems in 
each catalog such that GW networks will detect a sizable number 
of systems; the particular values of 4 x lO'^ NS-NS or 3 x 10^ NS- 
5 Mq BH systems are specified somewhat arbitrary. 



1983). In addition, in 



the actual Universe, we expect BHs in NS-BH systems to 
have considerable spin, and tidal disruption may occur 
for a NS orbiting prograde around a highly spinning 10 
M0 BH dFoucart et al.|201lD . We choose NS-NS and NS- 
BH systems with binary separations of 1.0 x 10~^Rq and 
1.4 X 10~^Rq respectively to ensure that they will merge 
within the system's characteristic gravitational radiation 
timescale (Shapiro & Teukolsky 1983). 

Regarding geometric source parameters, for each of 
our two catalogs, we distribute either 4 x 10^ NS-NS 
binaries or 3 x 10^ NS-BH binaries out to z = 0.5 
2.82 Gpc as suming a ACDM Universe given in Ko- 



matsu et al^2009 ) . Each binary in a catalog is associated 
with a random orientation such that p{cosi) ex const 
with cos L G [—1,1], and a random sky position such that 
p{cos6) (X const with cos 6 G [—1,1] and p{(j)) cx const 
such that G [0, 27r]. 

For those NS binary merger events with distances 
< 200 Mpc, we assume that the spatial distribution of NS 
binaries traces host galaxy light. We use a "Census of the 
Local Universe" (CLU) with information compiled from 
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particular galaxy is weighted by the B-band luminosity 
of that galaxy in CLU. The size of the galaxy is assumed 
to be three times the size given by the surface brightness 
contour at apparent magnitude 25 arcsec"^ in CLU. B- 
band luminosity incompleteness is taken into account by 
dividing the catalog into 10 Mpc bins and choosing ran- 
dom positions for galaxies which represent the missing 
luminosity. 

Finally, for those binaries located with distances > 200 
Mpc, we assume that the NS binary merger distribution 
has a constant comoving vol ume density in a ACDM Uni- 
verse (Komats u et al.||2009[ ). 

In summary, we construct two catalogs of 4 x 10^ NS- 
NS and 3 x 10"^ NS-5 Mq BH binary populations, where 
each binary is described by its set of physical and geomet- 
ric source parameters: {mi , m2 , , cos ^, cos 6>, cos (j)} . 

3. GW DETECTABILITY AND PARAMETER ESTIMATION 
- METHOD 

With the geometric and physical source parameters in 
hand for each binary, we can simulate the predicted GW 
strain emitted for every inspiraling NS binary in our two 
catalogs defined above. With knowledge of anticipated 
GW interferometric noise curves and by assuming ideal- 
ized Gaussian instrument noise, we can simulate the pre- 
dicted GW data-stream measured at a particular GW 
interferometer. Therefore, by matched filtering a GW 
detector output with a theoretically-predicted GW wave- 
form, measurements of GWs will allow us: i) to detect 
NS-NS and NS-BH binary inspiral and mergers, and ii) 
to extract the physical and geometric properties of the 
source. We first review our understanding of GW wave- 
forms, and then introduce the particular GW waveform 
that we use. Second, we outline the principles of matched 
filtering used when detecting and estimating parameters 
of the GW source. Finally, we describe the different 
GW networks considered and the three triggering cri- 
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teria used to construct different GW-detected NS binary 
merger populations. Further details can be found in NIO 
and Nil. 

3.1. GW waveform 

Turning to models of GW emission and dynamics, we 
view the GW waveform for merging compact binaries in 
terms of three phases; the inspiral, merger and ringdown. 
The inspiral phase, describing the loss of energy and an- 
gular momentum of the binary due to GWs, can be mod- 
eled accurately using the post-Newtonian (PN) approxi- 
mation in General Relativity. The PN approximation is 
an expansion in ~ jc? ^ where v is the characteristic or- 
bital speed for gravitationally-bound systems. The state- 
of-the-art accuracy for non-spinning inspiraling binaries 
is 3.5PN, corres ponding to an o rder of O (1/c^) in a PN 
expansion (e.g., |Blanchetl[2QQ6| . At 3.5PN, NSs and/or 
BHs are modeled using the 'point '-particle ('(^'-function) 
description, with finite-size effects being formally negligi- 
ble upto 5PN. However, several orbits prior to the merger 
of the two bodies, the weak-field PN approximation is no 
longer valid, and we require computationally-expensive, 
numerical simulations that model the merger pha se by 
directly solving Einstein Field Equations (see e.g. |Pre-| 
torius||2005|. After the two bodies have merged into a 
final single BH, perturbation techniques of a Kerr BH 
describe the ringdown. 

We make two important assumptions for the GW wave- 
form used in our work. Firstly, we assume that only the 
inspiral phase models the GW signal in this work. This 
is because the inspiral phase contributes to the major- 
ity of the signal accumulated in the frequency-band of 
advanced interferometers fo r NS-NS and several stellar- 
mass NS-BH systems (see Flanagan & Hughes 11998). 
Typically, the inspiral phase of NS binaries lasts &om a 
few to tens of minutes in the interferometer's frequency- 
band (Cutler et al.| |1993[ ). Secondly, we neglect the spin 
of NSs and BHs in our analysis. We expect NSs to have 
small spins in NS-NS and NS-BH systems. In contrast, 
BHs in NS-BH systems should have moderate to large 
spins. In the case of NS-10 Mq BH (with high spin) 
binaries, we expect spin precessional effects to increase 
the dimensionality of the parameter space and to modu- 
late the GW waveform significantly, which in some cases 
can improve GW sky localiz ation (van der Sluys et al. 



2008 Raymond et al. 2009). Consequently, neglecting 



the rnerger phase and assuming non-spinning binary sys 
tems will introduce systematic errors when estimating 
source parameters. However, as most GW detections 
will be at threshold, we estimate that statistical errors 
will dominate over GW waveform systematic errors for 
the majority of NS binary inspirals. 

The GW inspiral encodes a combination of the NS bi- 
nary's physical and geometric properties such as its red- 
shifted masses, its luminosity distance, its orientation, 
its source position, as well as the time and phase of co- 
alescence. Specifically, the predicted GW waveform at 
a particular detector a comprises the linear sum of the 
two GW polarizations /i+ and /ix weighted by the two 
antenna functions and Fx , and is given by: 



where D'^^ is the detector's response tensor, and de- 
notes the detector's position in spherical polar coor- 
dinates from the Earth's center. The scalar product 
n • denotes the time-of-fiight of the signal from the 
source to the GW interferometer. In our work, we use 
a GW waveform in the frequency domain, h-^{f) and 
hx{f)^ where the stationary- phase approximat ion as- 
sumes that {df/dt)/f <C / dProz et al l[T999| . The 



GW waveform used is accurate up to 3rN order in 
its phase (this improves the accuracy of the waveform 
used in NIO and Nil), and Newtonian in its ampli- 
tude. The 3.5PN GW phase depends on physical source 
parameters such as the redshifted chirp mass A4z = 

{l^z)Mc = (1 + ^)mi/^m2^V(^i + ^2)^/^ ^ is 
the binary's redshift (henceforth, the notation A4c refers 
to the physical non-redshifted chirp mass), and jHz = 
(1 + 2:) mim2 / {mi +7712) is the redshifted reduced mass. 
The 3.5PN phase also includes tc and <I>c, which are in- 
tegration constants and define the time and phase of 
coalescence respectively. In contrast, the Newtonian- 
order GW amplitude is a function of the GW frequency 
derivative / or so-called 'chirp' (itself a function which 
depends at Newtonian order on the A4z and at higher 
IPN order on the /i^), and of the geometric parame- 
ters such as the binary's cos^, Dl and its source po- 
sition (cos6>,0). Physical source parameters that ap- 
pear in the GW phase and / can be determined to a 
high accuracy for NS binaries because from thousands 
to tens of thousands of inspiral GW cycles could sweep 
up in the frequency band of the advanced GW detec- 
tors. However, only weak constraints on geometric source 
parameters, such as {cos Dl), are possible because of 
strong degeneracies that exist between parameters ap- 
pearing in the GW amplitude. This is the case for the 
majority of threshold-detected GW events (see NIO for 
further discussion). However, for (cos6>,0), differences 
in time-of-fiight among detectors in the network domi- 
nate over GW amplitude effects when reconstructing the 
event's sky position. We terminate our inspiral wave- 
form abruptly at the ISCO, /isco = (OV^ttM^)"^, where 
Mz = (1 + z){mi +777.2) is the redshifted total mass of 
the system. Such an abrupt cut-off of the GW waveform 
should have little impact on matched filtering for NS- 
NS binaries where /isco occurs at high frequencies with 
poor detector sensitivity (see Figure 1). In contrast, for 
higher total mass NS-BH binaries, /isco lies within the 
frequency band with high detector sensitivity and such a 
unphysical cut-off will introduce non-negligible system- 
atic errors. Finally, we ass ume that calibration measure- 
ment errors are negligible (Lindblom 2009 Vitale et al. 
20121). 
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3.2. GW Parameter Estimation 

Turning to GW parameter estimation, we summarize 
MCMC methods discussed in NIO and Nil. Our cen- 
tral quantity of interest is the posterior density function 
(PDF) of the distribution of inferred source parameters, 
denoted by the vector of parameters 6, following a GW 
measurement. The PN inspiral waveform used in our 
work depends on the vector 0, which comprises the pa- 
ranieters {M.z^ l^z, Dl^ cos 0, (j), cos l, iJj, t c , ^c}- Follow- 
ing 



Finn 



(n992) and [Cutler fc Flanagan' (199?) , 



we con- 



(1) sider a data-stream s{t) measured at a detector a that 
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comprises the instrument noise n{t) and a GW signal 
h{t,6)^ where 6 describes the source's "true" parame- 
ters, i.e., s{t) = n{t) + h{t^ 0). We assume that the noise 
at each detector has ideahzed Gaussian statistics. For a 
network of detectors, the PDF of the parameters given 
some set of observed datastreams s is: 



p(0|s)ocp(O)(0)£TOT(s|e), 



(2) 



where >Ctot(s | 0) is the total likelihood function and 
p^^\6) is the prior PDF that describes our prior knowl- 
edge of the signal's parameter distribution 0. The like- 
lihood function measures the relative conditional prob- 
ability of observing a particular datastream s given h 
and n. By assuming that the noise is independent at 
each interferometer, the total likelihood function £tot 
is equivalent to the product of the individual likelihoods 
at e ach detecto r. The likelihood Ca for detector a is given 
by ( [Ftol|[l992| : 



Ca {s I 6) = e 



= .-(^a(^)- 



\ha{e)- 



)/2 



(3) 



The notation {g\h) describes the noise- weighted cross- 
correlation of g(t) with h{t) in the vector space and is 
defined as 



{9\h) 



Jo 



df 



9*U)hU) + ~9U)h*{f) 



(4) 



Snif) 

where Sn{f) denotes the instrument's power spectral 



density. We discuss the form of Sn{f) in ^3.3 
Using Eqn. ([T]), for an ensemble of detector noise real- 



izations, the expected SNR at detector a is given by: 
S_ 

/ a, exp 

5" c 2 /GM ^^/^ 



96 Dl 7r2/3 V c3 
[F,%(l + cos^)^+4F,%(cos^)]^/^ 

/isco 



/ic 



df 



1/2 



(5) 



where /low is the instrument's low- frequency cut-off. Av- 
eraging over all possible binary sky positions and orien- 
tations, the expected sky -and-inc lination- averaged SNR 
is given by (see , Dalai et al.||2QQ6 ): 



sky- 



mc— ave 




For a binary that is directly face-on to an observer 



(cos^ = ±1), its SNRexp is a factor of v5/2 
1.58 greater than its inclination- averaged counterpart 
SNRinc-ave cit the samc position and distance, and a fac- 
tor of Y^5/4 1.12 greater than its s ky- and- inclinat ion - 
averaged SNRgky-inc-ave (see N 10 and Dalai et al.|2QQ6 ). 
Throughout the rest of the paper, we use Eqns. ([5]) 



and (|6| for GW-detectability (SNR-based) scaling argu- 
ments. We define the expected network SNR as the root- 
sum-square of the expected individual detector SNRs. 

To infer the geometric source parameters for 
each NS binary considered, in particular the subset 
{Dl^ cos cos 0^(j)), we map out the full posterior PDF 
of all source parameters using MCMC, given an observed 
data stream at a detector. The Metropolis-Hastings 
MCMC algorithm used is based on a generic v ersion of 
CosmoMC, described in [Lewis fc Bridle] ( [20021 ) . We as- 
sume prior distributions in all source parameters to be 
fiat over the region of sample space where the binary 
is detectable at an expected network SNR = 3.5. For 
each MCMC simulation used on a single NS binary in- 
spiral, we derive marginalized parameter measures and 
rms errors over {Dl, cos cos 6^ (j)) at 68%, 95%, and 99% 
confidence regions (henceforth, denoted c.r.). 

3.3. GW networks 

We consider GW networks consisting of combinations 
of LIGO (which comprises LIGO Hanford and LIGO Liv- 
ingston), Virgo, LIGO India, and KAGRA. In the rest 
of the paper, we use the following notation to describe 
different GW networks with n detectors: 

- Nets or network 3 is LIGO Hanford, LIGO Liv- 
ingston and Virgo, 

- Net4I or network 41 is LIGO Hanford, LIGO Liv- 
ingston, Virgo and LIGO India, 

- Net4K or network 4K is LIGO Hanford, LIGO Liv- 
ingston, Virgo and KAGRA, and 

- Nets or network 5 is LIGO Hanford, LIGO Liv- 
ingston, Virgo, LIGO India and KAGRA. 

Apart from LIGO India, the detector's positions (as 
measured from the Earth's center) used in this work are 
given in Table 1 of NIO. For LIGO India's position and 
orientation, we use East Longitude A = 76.7, North lati- 
tude = 14.3, Orientation = 0, x-arm tilt Vtx = ^ and 
?/-arm tilt Vty = {). 

For simplicity, we assume that the noise sensitivity 
curve for each detector is represented by the anticipated 
broadband-tuned sensitivity curve for a single advanced 
LIGO detector, shown in Fig. [l] We impose a low- 
frequency cut off at 10 Hz and frequencies below 10 Hz 
are not included in our analysis. In practice, LIGO, 
Virgo and KAGRA will have different noise sensitivities 
in different frequency bands because of variations in each 
instrument's design. We also consider the anticipated 
sensitivity curve for LIGO Livingston and LIGO Han- 
ford interferometers using optically squeezed light (|Miao| 
[2012|. To compare different GW network abilities, for 
each NS-NS or NS-BH binary, we assign a unique noise 
realization to each GW detector, which we keep constant 
when adding and subtracting detectors to a network. In 
addition, we assume that each GW interferometer oper- 



ates at an idealized 100% of the time (see Schutz||2QlT 
for estimates of different instruments' duty cycles 
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Fig. 1. — Anticipated noise curves for Advanced LIGO (solid 
black line; Harry Sz the LIGO Scientific Collaboration 2010) and 
Adva nced LiGU with optical squeezed light (dashed red line; Miao 
|2012| . In this paper, we assume a low frequency cut-off of lU Hz. 
i 'iie features at 10 and a few hundred Hz are various thermal noise 
resonant modes of mirror suspensions. The modes at a few hundred 
Hz are suspension fiber resonances. 



3.4. GW triggering criteria 

For each binary in the two catalogs, we implement 
three GW triggering scenarios that use: 

• a coincident trigger criterion (denoted "a"): we 
select the binary if at least two GW detectors 
each have an SNR > 6 and if the expected network 
SNR>12. The estimated number of GW tem- 
plates required and desired false-alarm rates deter- 
min es the ch oice in particular SNR threshold val- 
ues (|Owen|1996| and see e.g., Nil and ,Abadie et al., 
2QT2i| r 



• a coherent trigger criterion (denoted "6"): we select 
the binary if its GW expected network SN R> 8.5 
(see e.g.. Nil and Harry fc Fairhurst]|2Qll ). 



• an EM-precursor coherent trigger criterion: we 
select the binary if its GW expected network 
SNR > 7.5. Henceforth, to avoid confusion with 
the coherent trigger above, we refer to this case as 
an EM-precursor trigger criterion. As discussed in 
NIO, we choose to lower the network threshold in 
the presence of an already observed EM counter- 
part because prior knowledge of the merger time 
and sky position reduces the number of searched 
GW templates. Models of EM-precursor emission 
to NS binary mergers include resonant shattering of 
NS crusts o bservable in 7- and X-rays (e.g. |Tsang| 



et al. 



2012 ) and a coherent burst of radio emis- 

sion produced by magnetic ally d ominated outflows 
dPshirkov fcTostnov 2Q1Q| [PirQ|[2Q12 ). 



For all scenarios, the interferometers may have signif- 
icantly non-Gaussian noise, necessitating an increase in 
the SNR threshold to take into account complex detector 
statistics. 



4. GW-DETECTED BINARY POPULATIONS - RESULTS 

This section presents the relative fraction, rates and 
distributions of geometric parameters of NS-NS or NS- 
BH inspirals detected by different GW networks of in- 
terferometers. We deflne the term relative fraction to be 
the ratio of GW-detected binary mergers out of the total 
number of binaries in each catalog. 



4.1. Relative fractions of GW-detected events 

For each NS binary merger in our two catalogs, we 
compute and compare GW SNRs to a deflned thre shold 
SNR at a particular detector or at a network (§ 3.4). For 
different progenitors, GW networks, and t rigger mg cri- 



teria, we thus obtain fractions of those mergers that are 
detectable by GWs out of the catalog's total number of 
systems. Table[l]shows relative fractions of GW-detected 
merger samples. 

Table [T] also indicates the relative fraction of GW- 
detected mergers that have their orbital angular momen- 
tum vectors oriented towards the Earth such that they 
could show collimated 7- and X-ray emission. Observa- 
tions exist for tw o SGRBs indicating a half- jet opening 
an^le 0-^ of - 7° ([Burrows et"ar '2QQ6! 'Sode rberg et al.; 
2QQ6D and - 3 - 8^ ( [l^ong eT al. 2U12). A handful of 
other SGRBs exhibit upper and lower bound jet-break 
measurements, discussed in §5.1| In this work, we deflne 
beamed binaries to be those binaries whose orbital angu- 
lar momentum vector lies within a relatively stringent Oj 
of 6°. 

Table [T] illustrates several trends between different 
sampl es of GW-detected NS binary mergers. Schutzj 
( 2Q11[ ) provides powerful analytically-derived expressions 
that show good agreement with our explicit results. 
First, the fraction /beamed of beamed NS binary merg- 
ers seen in GWs, observable from all possible inclination 
angles ^, is less t han those binari es with isotropic orienta- 
tion (see |Schutz|2Qll| and Metzger fc Berger'2Q12 ). When 
^ <C 6>j, /beamed ^ 1 — COS 6>j ^ 6'^ /2 for Small ^j; in our 
case, the empirically-derived range /beamed ~ 0-5 — 1% 
agrees well with its theoretical value of 0.7%. Second, 
the relative fraction of GW-detected NS-5M0 BH merg- 
ers is greater by a factor from four to flve than GW- 
detected NS-NS mergers. This follows from Eqn. (Ig]), 

5 /6 

where SNR scales as Ale and detectable volume thus 
scales as mI^^^. For our NS-NS and NS-5M0 BH in- 
spirals, values for Aic are I.2IM0 and 2.22M0 respec- 
tively. Third, the fractions of GW-detected events where 
the two LIGO interferometers use squeezed light are a 
factor from ~ 9 to 10 greater than for those networks 
where no optical squeezing is implemented. Illustrated 
in Figure 1, the optically-squeezed advanced LIGO noise 
curve is a factor from two to three more sensitive than 
the standard analog's curve. From Eqn. ([6|, such an im- 
provement in instrument sensitivity translates to an im- 
provement by a factor of ~ 2^ — 3^ in detectable volume 
(because SNR is inversely proportional to Dl). Fourth, 
increasing the number of GW detectors in a network from 
3 to 5 increases the number of GW-detected mergers by 
a factor of ~ 2 or less (see NIO). Shown in Eqn. the 
network SNRexp scales as ^ y^, where n is the number 
of detectors. 
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4.2. Estimated relative rates of GW- detected inspirals 

To convert relative GW-detected fractions into relative 
GW-detected rate predictions^ we first require estimates 
of the astrophysical NS binary merger rate, independent 
of GW detection. In the case of NS-NS binaries, dif- 
ferent astrophysical merger rates are derived either by 
extrapolating the distribution of observed Galactic bi- 
nary pulsars or from population synthesis results (see 
e.g. , [Phinney'' 1991 and references in Abadie et al. 2010) . 
The rates range from 0.01 to 10 Mpc~^ Mvr~^, with 1 
Mpc~^ My r~^ being the mean of the rate's PDF (Abadie 
et al.||2010 ). In contrast, because we have yet to observe 
a NS-BH system, all estimates of NS-BH merger rates are 
based entirely on theoretical population synthesis results. 
The rates range from 6 xlO~^ to 1 Mpc~^ Mvr~^ , where 
0.03 Mpc~^ Myr~^ is defined as a realistic rate in Abadie 



et al. (2010). Therefore, our estimates for GW-detected 



merger rates rely on theoretical predictions of NS binary 
merger rates that span three orders of magnitude. 

Here, we estimate 7^ns-x, NS binary merger rates de- 
tected by networks of GW interferometers, using: 



A/ks-x X /ns-x X V X 



1 



(7) 



where the subscript X denotes a NS or BH, and A/ns-ns 
and A/ns-bh are the astrophysical NS-NS and NS-5 Mq 
BH merger rates in Mpc~^ yi"~^ respectively. V is the to- 
tal volume that the catalogs encompass (in our case, this 
corresponds to~4/3x7rx (2.82 Gpc)^) and /ns-x are 
the relative fractions of GW-detected NS binary merg- 
ers. Table 1 gives values of /ns-x for different progen- 
itors, GW networks and triggering criteria. The factor 
k ~ SVS applies to all netw orks with any number of 
detectors. It incorporates the | Abadie et al | pOlO) cor- 
rection for the GW interferometers' non- stationary and 
non- Gaussian noise, applied in order to achieve required 
false-alarm rates. 

Figure |2] shows relative rates of GW-detected NS bi- 
nary mergers which have either isotropic or beamed emis- 
sion. We use the m ean and /or realistic r ate of NS binary 
mergers quoted in Abadie et al. (2010). Given the few 
orders-of-magnitude uncertamty, we use NS-10 M© BH 
merger rates as representative for the merger rates of NS- 
SMq BH systems used in this work. For NS-NS mergers, 
we use a value of ^ 18 270 for the prefactor [A/ns-x ^] 
m Eqn. 0. For NS-BH mergers, we use a value of ~ 550 
for the prefactor [A/ns-bh ^j- Let us now discuss sev- 
eral features of Figure [2] 

First, out of all possible GW networks and triggering 
schemas, we present relative NS-NS and NS-BH binary 
merger rates detected by GWs under three scenarios: 

i) Nets a: coincident-triggered network 3 

ii) Net5b: coherent-triggered network 5 

iii) Net5b with optically-squeezed LIGO. 

Due to the high SNR threshold required at two detec- 
tors or more, NetSa detects the fewest number of NS- 
NS or NS-BH mergers. It provides a lower hound on 
the number of detected GW events, indicative of how 
the early years of GW measurements might unfold. In 
contrast. Net 5b detects the largest number of mergers of 



NS-NS or NS-BH mergers, because the coherent-network 
SNR scales as ^/n^ where n = 5 is the maximum num- 
ber of detectors. It hence provides an upper bound on 
the number of GW-detected mergers, suggestive of how 
a GW network might operate after the first several years 
of GW measurements. The third scenario, envisioned 
later in the timeline of the development of the GW 
network, provides a highly optimistic bound for NS-NS 
and NS-BH mergers detected using Net5b with optically- 
squeezed LIGO; we choose to investigate this scenario in 
a future study. Summarized in Table 2, we use NetSa 
and Net5b to indicate representative bounds for the per- 
formance between different GW networks and triggering 
schema. Illustrated by Figure |2j our GW-detected NS 
binary merger rat e estimates show go od agreement with 
earlier works (e.g., Abadie et al. 12 010]). Differences occur 
because of different triggering criteria invoked and SNR 
thresholds used. 

Second, the dark grey shaded regions in Figure [2] de- 
note those binaries that have their orbital angular mo- 
mentum vector lying within a relatively stringent 0^ < 6° 
and the light grey shaded regions denote those binaries 
whose orbital angular momentum have 0^ > 6°. From 
one to several NS binary mergers per year could have 
Oj < 6° and may exhibit 7-ray collimation associated 
with SGRBs. Our beamed NS-NS binary merger rates 
are consistent with SGRB rates of 10 Gpc~^ ^^7^^ 



cussed in l M etzger & Berger ( 20 12|), [Coward et al.| ([20121) , 
'Enrico Petf illo & Dietz (2012) and Chen & Holz'(2012). 



I'hird, we emphasize that lower and upper bounds to 
the rate estimates differ, for instance, for NS-NS merg- 
ers by two orders of magnitude below and an order of 
magnitude above from the mean values that we use. Re- 
sults presented in Figure |2] are instructive in that they 
illustrate relative GW detectability rates between differ- 
ent GW networks and triggering criteria, but the values 
given here should be used with caution. 



4.3. GW Malmquist effect in detected events' distance 
and inclination angle 

GW detection criterion sets implicit prior distribu- 
tions on geometric parameters of NS binary mergers 
(see NIO and Schutz 2011) . Defined here as the GW 
Malmquist effect, our G W detection criterion preferen- 
tially selects for more face-on (or equivalently beamed 
and more 'GW-luminous' binary inspirals); see Eqn. (|5|. 
The GW Malmquist bias is analogous to the standard 
Malmquist effect in observational astronomy, where in- 
trinsically brighter objects are detected further out. In 
GWs, beamed (cos^ ±1) binaries have higher SNRs 
and are i ntrins ically more luminous in GWs (Eqn. ([5])). 

Figure |3(a)| shows the 2-D distribution for parame- 
ters {Dl , cos l) of NS-NS mergers detected in GWs with 
NetSa. Important for EM follow-up and for coincident 
EM and GW observations, we remark on noteworthy fea- 
tures of the distribution. Figure [3(a) illustrates the GW 
Malmquist bias towards detection of beamed binaries, 
with cos^ ^ ±1. The distribution exhibits a character- 
istic V-shape which is consistent with the analytically- 
derived PDF of detected values in l given in Eqn. (28) 
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TABLE 1 

Relative fractions x10~^ of NS-NS and NS-5 Mq BH Mergers with collimated (denoted 'B') and isotropic 
(denoted T) emission Detectable in GWs Using Three Different Selection Criteria with Four GW Networks. The 
notation 'os' represents optical squeezing in the ligo interferometers. the range given represents the 1-a 

statistical error of our simulation. 



GW Network 




Nets 


Net4I 


Net4K 


Nets 






B|I 


B|I 


B|I 


B|I 


Coincident "a" 


NS-NS 
NS-5M0BH 


0.3 di 0.3|11 di 2 
0.7 di 0.5|50 ±4 


0.3 lb 0.3|17 ± 2 
1.0 ± 0.6|79 ± 5 


0.3 di 0.3|17 di2 
1.3 di 0.7|77 d: 4 


0.3 ± 0.3|23 ± 2 
2.3 ± 0.9|104 ± 6 


Coherent "b" 


NS-NS 
NS-SMqBH 


0.8 ± 0.4|36 ± 3 
2.3 ± 0.9|170 ± 7 


0.8 ± 0.4|57 ± 4 
3.7 ± 1.1|251 ± 9 


0.8 ± 0.4 |59 ± 4 
4.0 ± 1.2|243± 9 


1.5 ± 0.6|78 ± 4 
4.7 ± 1.2|323 ± 10 


EM precursor 


NS-NS 
NS-5M0BH 


0.8 ± 0.4|54 ± 4 
3.0 ± 1.0|244 ± 9 


1.0 ± 0.5|80 ±4 
4.0 ± 1.2|350 ±11 


1.8 ± 0.6|81 ± 4 
4.3 ± 1.2|350 ±11 


2.0 ± 0.7|113 ±5 
6.3 ± 1.5|464 ± 12 


Coincident V + O.S. 


NS-NS 
NS-5M0BH 


2.0 ± 0.8|129 ± 7 
6.7 ± 1.5|490 ± 12 


2.0 ± 0.8|140 ± 7 
8.0 ± 1.6|535 ± 13 


2.0 ± 0.8|140 ± 7 
8.7 ± 1.7|534 ± 13 


2.3 ± 0.9|152 ± 7 
9.3 ± 1.7|579 ± 13 


Coherent + O.S. 


NS-NS 
NS-5M0BH 


4.7 ± 1.2|364 ± 11 
27.7di3.0|1517 ± 21 


4.7 d= 1.2|391 ± 11 
30.0 dz 3.2|1643 di 21 


4.7 ± 1.2|390 ± 11 
30.7 ± 3.2|1640 ± 21 


5.0 ± 1.3|418 ± 12 
33.0 ± 3.3|1777 ± 22 


EM precursor + O.S. 


NS-NS 
NS-SMqBH 


8.0 ± 1.6|517 ± 13 
33.0± 3.3|2062 ± 23 


9.0 ± 1.7|565 di 13 
36.0 zb 3.4|2248 di 24 


9.3 ± 1.8|557 ± 13 
37.0 ± 3.5|2246 ± 24 


10.0± 1.8|610 ± 14 
39.3 ± 3.6|2425 ± 25 



TABLE 2 

Representative GW Network Scenarios for Detectable Samples of NS-NS and NS-BH mergers. The notation 

MED. AND MAX. REFER TO THE MEDIAN AND MAXIMUM VALUES OF PARAMETER DISTRIBUTIONS. 



Feature Lower Bound Scenario Upper Bound Scenario 



Relative fractions & rates Coincident 3 detector: Net3a Coherent 5 detector: Net5b 

19 yr-i (NS-NS) \ 3 yr'^ (NS-BH) 138 yr'^ (NS-NS) | 17 yr'^ (NS-BH) 

Detectable distance Coincident 3 detector: NetSa Coherent 5 detector: Net5b 

220-400 Mpc (med-max; NS-NS) 390-750 Mpc (med-max; NS-NS) 

350-600 Mpc (med-max; NS-BH) 650-1250 Mpc (med-max; NS-BH) 



Sky Area Errors 



Coherent 3 detector: NetSb 
55-180 deg2 (med-max; NS-NS) 
50-170 deg2 (med-max; NS-BH) 



Coincident 5 detector: Net 5a 
7-120 deg2 (med-max; NS-NS) 
6-65 deg2 (med-max; NS-BH) 
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Fig. 2. — Relative rate of NS binary mergers detectable by dif- 
ferent GW networks and triggering criteria. The dark grey shaded 
regions denote those binaries that have their orbital angular mo- 
mentum vector lying within a half-jet opening angle of 6° and the 
light grey shaded regions denote those binaries whose orbital an- 
gular momentum have half-jet opening angle greater than 6°. The 
notation 3 det.^ refers to NetSa, a coincident-triggered network 
3, 5 det. A denotes Net5b, a coherent-triggered network 5, and 
5 det. AO represents a NetSb with optical squeezing in the two 
LIGO interferometers. We use the mean and /or realistic rate of 
NS binary mergers quoted in |Abadie et al.| ( |2010| ). 



of Schutz (2011). Unsurprisingly, we detect the majority 
of events at threshold and observe a paucity of close-in 
binaries, detected with distances less than 100 Mpc. The 
maximum distance for NS binary mergers detectable by 
Net3a is - 400 Mpc. We note that the closest SGRBs 
with known redshifts are 080905, 050709, and 050724 at 
z = 0.122 (- 560 Mpc), z = 0.161 76 Mpc), and 
z = 0.257 (- 1.28 Gpc); see 'Berger| poTol ) . Therefore, 
the maximum detectable distance range of NetSa does 
not include the distances of the three closest SG RBs ob - 
served so far. Discussed in Metzger & Bergerl (|2012l), 
the lack of SGRBs observed within a few hundred Mpc 
is consistent with the Swift satellite's observational bi- 
ases: only ~ 1/10 of the sky is surveyed at a particular 
epoch and only ~ 1/3 of SGRBs observed by Swift have 
redshifts. Finally, in Figure 3(b)[ we show the 2-D distri- 
bution of NS-NS mergers for the parameters {Dl^cosl) 
detected by Net 5b. We note that the maximum de- 
tectable distance increases by a factor of 1.5 compared 
to NetSa. Moreover, in the case of our NS-5 M© BH cat- 
alogs, the SNR and hence detectable distance depends 

on the chirp mass MI^^ (Eqn. Therefore, the max- 
imum detectable distance increases to above 1 Gpc in 



Figure |3(c) 



4.4. Cumulative distribution of GW distances and 
localization errors 

Critical for EM follow-up, we examine cumulative dis- 
tributions in GW distance and sky area errors for pop- 
ulations of NS-NS and NS-BH inspirals detected using 
different GW detector networks and triggering criteria. 
For representative populations of GW-detected NS bi- 
nary mergers, we randomly take 200 NS-NS and 200 
NS-BH inspirals from their maximum sample detected in 
GWs with Net5b. We assume standard advanced LIGO- 
like noise curves with no optically squeezed light. 
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Fig. 3. — The 2-D marginalized prior distribution in Dl and 
cos L for GW-detectable NS binary mergers. Each point represents 
a GW-detectable NS binary merger. The top panel shows NS-NS 
mergers detected by NetSa. The middle panel shows NS-NS merg- 
ers detected by NetSb. The bottom panel shows NS-BH mergers 
detected by NetSb. Redshifts are computed assuming cosmological 
parameters given in ^Komatsu et al., ( ,2009, ). 



Figures 4(a) and 4(b) show the specific distribution 
in luminosity distance (in Mpc) of NS-NS and NS-BH 
mergers detected by different GW networks and trigger- 
ing schema. We use the term specific because we nor- 
malize the cumulative distribution to the sample of NS 
binary mergers that a particular triggere d G W network 
can detect. Shown in Table 2, similar to §4.1[ Net 3a and 
NetSb provide representative lower and upper bounds of 
the GW detectable distance. For NS-NS mergers, we 
find that median detectable distances are 180 Mpc and 
370 Mpc with Net3a and NetSb respectively. For NS-BH 
mergers, we find that the median detectable distances 
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are 240 Mpc and 660 Mpc with NetSa and Net5b respec- 
tively. 

Illustrated in Figures |4(a) and |4(b)[ two distinct dis- 
tributions for detectable distance exist depending on 
whether the GW trigger is coincident versus coherent. In 
contrast, we find that detectable distance ranges depend 
only weakly on t he number of detectors in a network. 
From Figure [4(b)| the detectable distance ranges for NS- 
BH mergers are approximately a facto r of two greater 
than for NS-NS mergers (see Eqn. ([6|, ^Oand |43|. 

In addition to distance, EM follow-up d etecta bilit y re- 



lies on sky area error ra nges; see Figure s 5 (a 



and ^ From Figures 5(a) and 5(b^ 



we 



and 5(b) 
find that a 

coherent-triggered network 3 provides the largest 
sky area errors for GW mergers. We refer to this scenario 
as NetSb and it represents our 'lowest-bound' on sky area 
errors (Table 2). On the other hand, a coincident- 
triggered network 5, denoted as NetSa^ provides the 
smallest sky area errors and represents our 'upper bound' 
for sky localizat ion (T abl e 2). 
From Figures [5(a)l and [5(b)| we find that 50% of NS- 



we 

.2 



NS mergers are detected to within 7 deg with Net5a and 
to within 60deg^ with Net3b. For NS-BH mergers 
find that 50% of events are detected to within 6 deg 
with Net5a and to within 55deg^ using Net3b. As ex- 
pected, similar distributions in sky area error exists be- 
tween GW-detected NS-NS and NS-BH merger popula- 
tions because most events are detected at threshold SNR. 

We find elliptically-shaped sky errors for the major- 
ity of our examined NS binary mergers (see ^6.1| for an 
example); differing GW arrival times at each detector 
dominate sky area reconstruction rather than parameter 
degeneracies in the GW waveform's antenna functions 
(see Nil). In a handful of cases, we find multimodal 
peaks for especially weak SNR events because of larger 
uncertainties in arrival times at detectors. In addition, 
we find that sources located in (or close-by to) the degen- 
erate Ne t3 plane have relatively poor angul ar res olution 
^irhurst||2010| 



(see a lso 



et al. 



20 



Nil, Veitch 



Wen fc Chen|[2QTQ[ 

An improvement by a factor of two in the 

normalized cumulative sky error is seen with network 41 
(LIGO+Virgo with LIGO India) compared to network 
4K (LIGO+Virgo with KAGRA) only in the case of using 
a coincident trigger. Given that LIGO India is located 
further away from the degenerate LIGO- Virgo plane than 
KAGRA, such a facto r of two impr ovement in sky area 
error is expected (e.g., SchutzpOll ). 



4.5. GW volume estimates 

Measurements by GW networks provide us with dis- 
tance and sky area errors. With both values in hand, we 
can construct GW volumes, which aid in identifying the 
EM counterparts of NS binary mergers (as ^6.2 describes 
in detail). 

As a first attempt, we introduce and define below the 
term low-latency GW volumes. Such volumes in prin- 

^ In contrast, Figures 3 and 4 in Nil show sky error distribu- 
tions for subsets of NS binary mergers detected by different GW 
networks that are normahzed to the fuh detected sample by net- 
work 5. In this work, instead of emphasizing the reduced number 
of detections, we particularly wish to answer what % of NS-binary 
mergers are detected by a known triggered network to a certain 
sky area error. 
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Fig. 4. — Cumulative luminosity distance (in Mpc) distribution of 
detected sample of NS-NS (top panel) and NS-BH (bottom panel) 
mergers normalized to each specific network and trigger criterion. 
The blue lines denote those NS mergers detected using a coincident- 
trigger criterion; the red represent those events detected triggering 
using the GW network coherently. Solid lines represent GW net- 
work 3, dotted lines denote GW network 41, dashed are GW net- 
work 4K, and dash-dot are GW network 5. 

ciple can be computed within a few to tens of minutes 
of a GW detection (and do not rely on the fuh MCMC 
machinery used in this work) and are hence critical for 
EM follow-up. In their final science run before halting 
for upgrades to their advanced versions, LIGO and Virgo 
sent triggers to EM telescopes within ~ 30 minutes of a 
possibl e GW signal bei ng detected 



, „ (LIGO Scienti fic Col- 

labo ration et al.|2012a[ Abadie et al. 2012b; Evans et al. 



|2012| ; most of this time was spent for human- limited 
verification checks at each detector site. In the era of 
advanced detectors, efforts are underway t o reduce the 
latency timescale to less t han ten minutes (Singer et al. 
2012[ ICannon et al.|[2QT2l ). 



Tn this work, we compute low-latency GW volumes 
by using only marginalized 2-D sky area errors and 
marginalized 1-D distance measures (all at 95% c.r.). As 
we now discuss, although in this work we derive sky area 
and distance errors by marginalizing the full 9-D PDF, 
we could instead have used ~ minutes timescale compu- 
tations of approximate sky area errors and distance mea- 



Identifying the EM Counterparts of GW Mergers 



11 



c 

DD 

"D 
CD 
■*—> 

o 

CD 
■*—> 
CD 

Q 



o 

CD 
Q. 
CO 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



Coincident 



Coherent 



10 100 
Area (sq deg) 

(a) NS-NS binary mergers 




10 100 
Area (sq deg) 

(b) NS-BH binary mergers 



Fig. 5. — Cumulative sky error (in deg^ at 95 % c.r.) distributions 
of detected sample of NS-NS mergers normalized to each specific 
network and trigger criterion. The blue lines denote those NS-NS 
(top panel) and NS-BH (bottom panel) mergers detected using a 
coincident-trigger criterion; the red represent those events detected 
triggering using the GW network coherently. Solid lines represent 
GW network 3, dotted lines denote GW network 41, dashed are 
GW network 4K, and dash-dot are GW network 5. 

sures. Regarding sky localization errors for the majority 
of G W- detected mergers , analytically-derived formulae 



(e.g. Wen fc Chen 2010 



the seconds timescales, al 



Fairhurst 2011), computed on 
. ow for sky reconstruction esti- 
mates that are in good agreement with explicitly-derived 
2D sky errors presented in Nil. This is because differ- 
ent GW arrival times at each detector dominate over 
amplitude corrections in the GW waveform. Regarding 
distance measures, Fisher matrix-based estimates allow 
for rough distance mea sures on a second timescale (see 



Ajith fc Bose|200"9 ). In practice, however, measured 



e.g. 

distances lor the majority of threshold events will have 
significantly larger errors (by a factor of a several) from 
their Fisher-matrix derived counterparts (NIC). This is 
because degeneracies between the sources' geometric pa- 
rameters that appear in the GW waveforms' amplitude 
inhibit measurement inference for low SNR events (see 
NIO for a detailed discussion). Therefore, a possible so- 
lution when estimating low-latency distance measures is 
to use their Fisher matrix-derived errors multiplied by a 



factor of three (see NIC an d [Del Pozzo||2012D . 

To compute low-latency GW volumes, we use upper, 
mean and lower distance measures; we define d^ and d/ to 
be the upper and lower 1-D marginalized distance values 
at 95% c.r.. We replace d^^ with the horizon or max- 
imum detectable distance of a coincident- or coherent- 
triggered GW network when d^^ > d^. We show absolute 
volumes in Mpc^ as a function of the mean distance with 
their upper and lower distance errors for NS-NS and NS- 
BH mer gers d ete cted b y Net3a and Net5b respectively 
(Figures 6(a) and |6(b) ). As expected, the GW measured 
upper and lower distance ranges are noticeably smaller 
for those NS-NS and NS-BH binaries with true distances 
less than 200 and 500 Mpc respectively. 



5. EM DETECTABILITY 

In this section, we first review characteristics of sug- 
gested EM counterparts to compact binary mergers. We 
then discuss EM detectability by upcoming or current 
optical and infrared telescopes. 

5.1. Predicted EM counterparts to NS binary mergers 

Short gamma ray bursts (SGRBs): The leading pro- 
genitor models for the majority of o bserved SGRBs are 
NS-NS and NS- BH mergers (see e.g. [E ichler et al.||1989" 

Pr^ 1 • -1 rkml Fivr j. i lli rkrirft r i i "r\ 



[Paczynski 1991] |Narayan et al.||1992p. The hypothesis 
nas been further supported by around forty SGRB ob- 
servations triggered by the Swift satellite and fo llowed- 
up by rapid multiwavelength observations (e.g. Berge 



et al 



et al 



2005 



2006 



Fox et al.|[2QQ5l [Hjorth et aLl|2QQ5 
• Berger||201lf 7"^" 



Bloom 

Theoretical models assume 

that accretion by a rotationally-supported disk onto a 
newly formed BH (or rapidly rotating NS) powers a rel- 
ativistically collimated outfiow, which res ults in the ob- 
serve d prompt gamma-ray emission 

_ 



. Rosswo g et al. 
'Rezzolla et al. j|201ip .' 

and energies, tne pror 



[e^ 



Ruffert et al. 



20031 IShibata fc l'amguclil[|:^0'n^ 



Due to their high Lorentz factors 
prompt 7-ray emission is assumed to be 
relativistically beamed with initial gamma-ray emission 
that lasts for < 2s (hence th e use of the name "short" 
when classifying SGRBs; see Nakar 2007 for a review). 
Jet-break observations in at least two SGRBs suggest col- 
limation of half opening angles of ~ 7° ( Soderberg et al.| 
2006[ [Burrows et ar| |2Q06) and 3-8° (Fong^etalJ[2^r 



r esp ectivel y. Upper and lower limits exist in^ lew other 
cases (e.g. [Fox et al.||2005| iGrupe et al ||2006[ ). As the 

relativistic beamed outtLow interacts with the surround- 
ing medium, we expect to observe afterglow signatures 
in the X-ray and optical occurring at longer timescales 
from minutes to days. Observations of EM afterglows 
suggest energies of E < 10^ ^ erg and circumburst densi- 



ties o f n < 0.1 cm ^ (Berger et al.|200'5 Soderberg et al 



20061). Afterglow model predictions as a lunction ol E 
and n are given in [van Eerten fc MacFadyen (2011); 
Eerten et ol] ( |2012| : 

K-process radioa ctivity transient s-kilonova\ Initially 
proposed by Li fc Paczynski (1998), rapid (r)-process 
radioactivity-powered transients are weak supernovae- 
like events. In this paper, we refer to such transients 
as 'kilonovae', so-called because their predicted peak lu- 
minosities are estimated t o be a fa ctor ~ 10^ greater 
than standard novae ( ^Kulkarni|2005| ,Metzger et al.^20TQ{ 
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Fig. 6. — Absolute measures in volume (Mpc"^) for a detected 
sample of NS-NS mergers observed by GW Net3a (top panel) and 
NS-BH mergers observed by GW Net 5b (bottom panel). Each filled 
point represents a detected NS binary merger at is mean luminosity 
distance. The horizontal error bars represent the upper and lower 
1-D marginalized distance values at 95% c.r. for each NS binary 
merger. The different col ors r epresent different case studies of NS- 
NS mergers examined in ^6.1| green is Case I (beamed binary), red 
is Case II (a nearby binary j, fight blue is Case III (a merger at low 
Galactic latitude), purple is Case IV (a merger at high Galactic 
latitude) and orange is Case V (a binary in a dense galaxy cluster 
environment). Specifically: green square is Case I, red circle is 
Case II, light blue triangle is Case III, purple diamond is Case IV, 
and orange star is Case V. 



[Roberts et al. 2011). A central premise of the model is 
tliat NS mergers produce ejecta from either dynamically- 
ejected tidal tails, or accretion disk outflows driven by 
early neutrino winds or late thermonuclear driven winds 
(e.g. 'Kul Wni|[2QQ5l [Metzger et aL| |2QQ8[ |2QQ9[ [Dessart 
|et al., 2QQf| r 



does not fall back onto the newly-formed BH or rapidly 
rotating NS. Numerical relativity and SPH simulations 
predict two orders of magnitude difference in the mass 
of the ejecta (0.001-0. IM©) and a fact or of a few in 



the ejecta's ve locity (0 . 1-0.3 c ); see e.^., Rosswog et al 
(1999); .Ra ntsiou e t al.| (|2008 |); i Foucart et al.| ( |201l[ ); 



ran et al 



(2012); East et al. (2012a|bp. Following tEe 



expansion of neutron-rich material from nuclear densi- 
ties, r-process nucleosynthesis produces heavier unsta- 
ble radioactive elements which subsequently beta-decay 
and fission back to stability on longer timescales. We 
expect the material to act as a heat source, and the 
subsequent emission to radiate isotropically. Based on 
highly uncertain opacities, light curves and color evo- 
lutions using radiative transfer models suggest that the 
emission peaks either in the optical or near infra-red. 
In the optical, the emission could peak with luminosi- 
ties of 10^^ — 10^^ ergs ~^ which decay o n the half to 
five day timescales ( [Metzger et al. [2010 ). Peak abso- 
lute magnitudes Mji range trom —14 to —17 mag and 
depend on the assumed ejecta mas s, velocity, opacity 
calcu lations and nuclear reactions (Metzger & BergerJ 
|201 2|). On the other hand, preliminary work estimate 
peak absolute magnitudes Mh = —15.5 mag in the near- 
infrared (~ 1.7/im) assuming an ejecta mass of O.OIMq 
at 0.1c with timescales varying from several to tens of 
days timescale (Kasen 2012). Efforts are currently un- 
derway to predict the spectroscopic Doppler-broadened 
signature of kilonovae. 

Radio counterparts: There are three predicted radio 
counterparts: i). we expect observable non-thermal radio 
emission from beanied ultra-relativistic ejec t a of SGRBs 
(|Berjg:er et al.l[2005[ [Soderberg et al.[|2006| [Chandra fc 
l ^mil||2012[ ), iij. we could observe a conerent radio burst 



emitted from a magnetically-d riven, relativist ic plasma 



outflow prior to the NS me rger ( Hansen & Lyutikov 2001^ 
Postnov|2Q10 ), and iiij. recent work suggests 



Pshirkov . 



he ejecta is gravitationally-unbound and 



incoherent radio signatures from blast waves produced 
by the interaction of quasi-spherical, sub- o r mildy- 
relatiyistic ejecta with the interstellar medium ([Nakar &| 
[Piran''2011 ). Nakar & Piran (2011) estimate tfiaJ radio 
nares may peak at 1.4 GHz emission for weeks out to red- 
shifts of 0.1 (~ 450 Mpc), and can be detectable at mil- 
liJansky levels. If the outflows are sub-relativistic, flares 
may be detectable on the years timescale at 150 Mpc at 
closer distances with current and near-future surveys. 

5.2. Differences between EM counterparts 

We highlight four features that distinguish the pro- 
posed EM counterparts and help define search strate- 
gies. First, the counterparts exhibit either beamed or 
isotropic emission. SGRBs have collimated jet emission 
and only accompany a very small fraction of NS-NS and 
NS-BH mergers (Table 1). On the other hand, kilono- 
vae and radio remnants have predicted isotropic emission 
and accompany all NS-NS and NS-BH mergers. Sec- 
ond, there is a wide disparity in timescales for fast and 
slow counterparts. SGRBs last for seconds (and their 
afterglows decay as a power law in time), kilonovae for 
hours to days and radio transients for months to years. 
Third, the rate of false positives is considerably differ- 
ent across EM wavelengths. The precise timing of the 
SGRB overcomes the challenge of poor sky localization 
of GW events. The quietness of the transient radio sky 
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is a boon to t he small number of spatially coincident 
false positives (Frail et al. 2012). The dynamic optical 
sky results in tens to hundreds of false positives that 
would be spatially and temporally coincident with GW 
events and search strategies are necessary to separate the 
wheat from the chaff. Fourth, discovery and follow-up of 
SGRBs is now a mature field. For the handful of mergers 
beamed towards us, we have rehearsed what needs to be 
done. On the other hand, off-axis and orphan SGRB af- 
terglows, kilonovae and radio transients are an uncharted 
territory. Both observational and theoretical progress is 
ongoing in leaps and bounds as we prepare for GW de- 
tectors to come online. Theoretical models continue to 
become more sophisticated with their predictions. Ob- 
servationally, synoptic surveys in the optical are already 
uncovering entirely new classes of fainter and rarer tran- 
sients. A suite of new radio facilities and radio transient 
searches are also coming online. For instance, wide-field 
low frequency (say, < 1 GHz) radio detectors (e.g. LO- 
FAR, MWA, EVLA) should be sensitive to pre-merger 
coherent emission and coincident timing can be used to 
connect them to GW detections. Relatively higher sensi- 
tivity and higher frequency radio detectors (e.g. EVLA, 
ASKAP, Apertif) are well-suited to searching for radio 
relic emission months to years after the GW detection. 

5.3. Detectability of EM counterparts 

Next, we discuss the detectability of isotropic EM 
counterparts by optical and infrared telescopes. As quan- 
tified in §4.4, the localization and distance horizon dis- 
tribution are dependent on the number of detectors and 
the threshold criterion in the GW network. We consider 
here the extreme cases for GW maximum detectable dis- 
tances (Net3a and Net5b), and for sky error areas (Net3b 
and Net5a). 

5.3.1. Optical Facilities 

Astronomers have a diverse arsenal of optical tele- 
scopes worldwide. We limit the discussion here to only 
telescopes with cameras larger than 1 deg^ (given the 
large localization areas) and apertures larger than Im 
(given the faintness of predicted counterpart). We con- 
sider current or scheduled-to-be operational telescopes. 
We divide telescopes into three categories: lm-3m class 
telescopes, 4m- 7m class telescopes, 8m- 10m class tele- 
scopes. Table |3] provides a summary of sensitivity and 
FoV of each telescope and camera system. 

Theoretical predictions of optical EM counterparts 
span orders of magnitudes in both predicted luminosity 
and predicted timescale (§ 5.1). To evaluate the relative 
merits of follow-up with dinerent telescope facilities and 
to begin to define a search strategy, we need to make 
a conservative assumption on the nature of the counter- 
part. Hence, for the discussion below, we first assume 
that the optical counterpart of a NS-NS merger will be 
brighter than Mr = —14 mag for at least two hours 
(later, we relax this assumption to —11 mag). Let us say 
a particular telescope takes three images at a separation 
of one hour, and each of these images has a 5-cr depth of 
Mr = —14 mag. Then, the counterpart will be discov- 
ered with high SNR (e.g., 12. 5-cr if the transient peaked 
at Mr = —15 mag) in the first image, at least at 5-cr in 



the second image and possibly below a 5-cr threshold in 
the third image. This is our "minimum" criterion for a 
secure detection. If the counterpart is either more lumi- 
nous or evolves at a slower rate, it will only improve the 
security of our detection. We require a minimum of two 
detections to securely distinguish the optical counterpart 
from moving objects in our solar system (asteroids) and 
artifacts. 

Telescope time is a zero sum game. A telescope with 
a given FoV of camera and a given aperture has to per- 
form a three-way tradeoff between depth, cadence (how 
frequently the same field is observed), and area covered. 
Here, we assume that each telescope takes at least three 
sets of images separated by one hour. In this one hour, to 
attain the Mr = —14 mag sensitivity for the most num- 
ber of events, the telescope will either integrate longer on 
a given field to see events further away or map a larger 
fraction of the localization area. If a telescope has a large 
aperture and small field camera, it will spend the one 
hour taking short exposures on a larger fraction of the 
localization area. If a telescope has a small aperture but 
a large field camera, it will spend the one hour stacking 
images to maximize integration time. It is precisely this 
choice that determines how many optical counterparts 
are detectable by a given telescope {modulo idealized ob- 
serving conditions as we shall discuss). 

We quantify the implications of this trade off on the 
number of detected NS-NS mergers in Figure 7(a) for a 
GW Net5b. 

For example, let us consider the r ole of CFHT in GW 
Net5b (green open squares in Figure 7(a)). In 100 s (60s 
exposure + 40s readout), CFHT can take a 0.9 deg^ im- 
age with a depth of 23 apparent mag. In one hour, 
CFHT can take 36 exposures, hence there are 36 pos- 
sibilities for the tradeoff. We discuss the first and last 
point on the curve of green squares. If CFHT spent the 
entire hour integrating on only one field, it would achieve 
a depth of 24.9 apparent mag and detect binaries with 
distances less than 615 Mpc (99%) but localization ar- 
eas less than 0.9 deg^ (2%). Instead, if CFHT spent the 
entire hour covering the large localization area and only 
spent one minute per field, it would achieve a poorer 
depth of 23 apparent mag and detect binaries only out 
to less than 250 Mpc (20%) but localization areas less 
than - 32deg2 (73%). 

Next, let us take the case of 8m-class Subaru's Hy- 
perSuprimeCam (HSC, blue open squares). By spending 
only five exposures on a given field, the depth of HSC 
can cover 100% of distances of detected NS binary merg- 
ers. However, it's smaller FoV camera limits the total 
area covered in one hour to 127deg2 i.e. 96% of mergers. 
LSST has the same depth but a larger camera, it can 
detect 100% of mergers. 

Let us next st udy t he case of Im-class ZTF (red filled 
circles in Figure 7(a) ) in a GW Net5b. In less than five 
pointings, with its superior 35 deg^ camera, ZTF can 
cover 100% of all localization areas. But its small aper- 
ture limits sensitivity to apparent 22.6 apparent mag or 
210 Mpc i.e. 10% of the mergers. 

Now, we consider the implications of the isotropic op- 
tical counterpart being much less luminous, for instance 
Mr = —11 mag (Figure [7(b) I and Table [4|. The percent- 
age of detectable counterparts goes down from 100% to 
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TABLE 3 
Optical Telescopes 



Telescope 


Aperture 


Field of View 


Exposure 


Overhead 


Sensitivity 


Reference 




(m) 


(deg^) 


(sec) 


(Readout) 


(5(7 mag in R band) 




^.'xiiT^ r*w"\r T^yj^ "n Qi rf^^Trl" Ti^j^r^i n'l"Ar 

/J W Iv/xv y X 1 diioldilj -L CtVylilLy 


1.2 


35 


60 


15 


20.6 


a 


La Silla Quest 


1.0 


9.4 (80%) 


60 


30 


20.5 


b 


Catalina Real-Time Transient Survey 


0.7 


8.0 


30 


18 


19 


c 


Palomar Transient Factory 


1.2 


7.1 


60 


40 


20.6 


d 


Pan-STARRS 1 


1.8 


7.0 


60 


3 


22.0 


e 


Skymapper 


1.35 


5.62 


110 


20 


21.5 


f 


CTIO-Dark Energy Camera 


4.0 


3.0 


50 


17 


23.7 


g 


WIYN-One Degree Imager 


3.5 


1.0 


60 


30 


23 


h 


CFHT-Megacam 


3.6 


0.9 


60 


40 


23 


i 


Large Synoptic Survey Telescope 


8.4 (6.7) 


9.6 


15 


2 


24.5 


J 


Subaru-HyperSuprimeCam 


8.2 


1.77 


30 


20 


24.5 


k 



^ iHadjiyska et al. |2Qll| 



^ Drake et aL 2 Qn9r 
Law et al. 2009 



see http: / / pan-starrs.ifa.hawaii.edu" 

^ see http:/ /rsaa. anu.edu.au/observatories/siding-spring-observatory/telescopes/skymapper/skymapper-instrument 



S Bernstei n et al.l20l3 

^ see http: //www. wiyn.org/ODI/Obsery^^wi^nodiove^i^ 



^ see http: / /www .cfht.haw aii.edu/Instruments/Imaging/Megacam/generalinformation.html| 
j IlSST Science Collaborations et aLhoOol , 

^ see http://www.naoj.org/Projects/HSC/index.html and http://www.naoj.org/cgi-bin/img_etc.cgi 



82% for LSST, from 96% to 42% for HSC and 97% to 
16% for DES. 

Assuming the most optimal strategy is chosen for each 
merger in this simulation, we can compute the fraction 
of detectable optical counterparts by each telescope (Ta- 
ble |4|. As exemplars, the smallest and largest FoV cam- 
era m each telescope aperture-class is chosen. Initially, 
when there is a GW three-detector network, binaries 
would be detected closer in and the localizations would 
be poorer. The smaller telescop es wit h larger FoVs will 
play an important role (Figure 
GW five-detector network, once 



a) ). In the era of a 
realization is improved 



and maximal detectable distance pushed fu rther back, 
the larger telescopes will be essential (Figure [8 (b)[ ). 

Finally, we consider the case of NS-BH mergers. Given 
that NS-BH mergers will on average be detected a factor 
of two further away, but have predicted optical coun- 
terparts 1.5 mag brighter, we get similar detectability 
fractions as NS-NS mergers (Table J4|. 

We emphasize that the detectable fractions presented 
in Table ffl are relative and subject to two caveats. First, 
there would be tiling inefficiency and edge effects due to 
the irregular shapes of GW localization and the rectan- 
gular/circular fields of view of the EM cameras. Second, 
all optical telescopes in this discussion are subject to cer- 
tain reality checks — they cannot observe too close to the 
sun or too close to the moon, if it is cloudy or raining, 
or if the target is in the quadrant of sky not accessible 
from a given location. Typically, these factors amount 
to^x|x| = |of the targets being visible at a given 
telescope on a given day respectively. 

We conclude that a network of telescopes at different 
longitudes, latitudes and mountain-tops would maximize 
the odds of follow-up. Hence, the numbers presented here 
should only be interpreted as illustrative of the relative 
detectability by different telescopes. 



5.3.2. Infrared Facilities 

Recent theoretical calculations of kilonovae opacities 
suggests that a significant fraction of the luminosity may 
be em itted in the redder bands beyond 1 /im ( |Kasen| 
|2012|). Unfortunately, our current suite of near-infrared 
facilities is not as wide-field as the optical with no camera 
larger than a square degree. 

Currently, the two widest field infrared facilities are the 
0.594 deg2 VIRCAM on the 4.1m VISTA telescope and 
the 0.19 deg^ WFCAM on the 3.8m UKIRT telescope. 
Fortunately, efforts are underway to build a 6.5m SASIR 
telescope with a 0.2-1 deg^ camera (SASIR,2012J . More- 
over, unlike VIRCAM and WFCAM, SASiR is expected 
to have a contiguous focal plane and simultaneously im- 
age in YJHK-bands. 

Efforts are also underway to build two wide -field in- 
frare d satellites — WFIRST ( |Green et al.pQT2| and Eu- 
clid dAmendola et al.||2012D . A set of three WFIRST 
fields is expected to cover 2 deg^ and eighteen minutes 
per field will give a 5-cr depth of H ~ 25.9 mag for imag- 
ing and H ~ 23.5 mag for low-resolution prism spec- 
troscopy. Euclid would need four fields to cover 2 deg^ 
and achieve depths of H ~ 25.6 mag and H ~ 21.4 mag 
in grism mode in the same time. Both missions a re con- 
strained to observe ~ 90° away from the Sun (Hirata 
eraI][2QT2| . 



We remmd the reader that the median-maximum lo- 
calization of Net 3b is 55-180 deg^ and Net 5a is 7-120 
deg^ (Table 2). Hence, the infrared follow-up will require 
an extremely large number of pointings to tile the area 
and will be limited to the best localized binaries. 

6. IDENTIFYING EM COUNTERPARTS 

Detection of candidate EM counterparts is only the 
first step. The most pressing question for EM telescopes 
looking at vast sky areas will be whether the transient 
objects are true GW emitters or false-positive signals 
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TABLE 4 

Relative percentages of isotropic, optical counterparts of GW-detected mergers detectable by different telescopes as a function of GW 

NETWORK, TRIGGERING CRITERION AND PEAK OPTICAL LUMINOSITY. INCORPORATION OF REALISTIC OBSERVING CONDITIONS (mOON, SUN, WEATHER, LATITUDE 

ETC.) REDUCES EFFICIENCY BY ~ 1/4 









PTF 


ZTF 


CFHT 


DES 


HSC 


LSST 


NS-NS merger 


& Mr < 


— 14 mag 














GW Net3a 






39 


44 


69 


97 


97 


100 


GW NetSa 






34 


41 


9S 


98 


98 


100 


GW Net3b 






18 


22 


34 


82 


79 


100 


GW NetSb 






10 


13 


61 


97 


96 


100 


NS-NS merger 


& Mr < 


— 11 mag 














GW Net3a 












19 


39 


86 


100 


GW Net5a 












18 


48 


91 


100 
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8 


16 


4S 


93 


GW Net5b 
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82 
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— IS.S mag 














GW Net3a 






66 


79 


79 


97 


97 


100 
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63 


72 


93 
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93 
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24 


39 


36 


78 


76 
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GW NetSb 






22 


28 


S6 


96 


94 


100 


NS-BH merger 
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— 12.S mag 














GW Net3a 






3 


3 


17 


79 


93 
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GW NetSa 






2 


2 


IS 


78 


98 


100 


GW NetSb 






1 


1 


S 


33 


47 


98 


GW NetSb 






1 


1 


S 


30 


S3 


88 



mimicking an EM counterpart. For instance, the optical 
sky is so dynamic that there wih be hundreds of fore- 
ground and background false positives associated with 
any detection. Foreground signals are, for example, M- 
dwarf flares, CVs and other stellar variables in the Milky 
Way. The foreground rates, therefore, depend strongly 
on the Galactic latitude and have a wide range of ampli- 
tudes and timescales. Background signals are supernovae 
(SNe) and AGNs at higher redshift than the GW de- 
tectable distance horizon for NS binary mergers of ~ 200 
Mpc - 1 Gpc. Thanks to systematic optical synoptic sur- 
veys, rate estimates of different classes of SNe occuring 
in a range of galaxy hosts now exist; for instance, core- 
coll apse SNe rate is 7.1 (±0 l)x 10"^ Mp c~^ yr"^ (see 



e.g., [Leaman et aL||2011| |Li et al.||2Qlla|bD . 

Hence, panchromatic follow- up (especially optical 
spectroscopy) is critical to unambiguously associate the 
counterpart with the GW signal. Given predicted opti- 
cal light curve evolutions, the timescale for spectroscopic 
follow-up should be within the hours to day timescale. 
There are a large number of telescopes in the 3-5m class 
range which can easily take low resolution spectra of 
transients brighter than 21 apparent mag. However, op- 
tical counterparts will likely be in regime where the tran- 
sient is fainter than 22 apparent mag and a > 6m-class 
telescope will be needed for spectroscopy. The list of 
such telescopes is rather small: the twin Magellan 6.5 
telescopes, the MMT 6.5m telescope, the twin Gemini 
8m telescopes, the four VLT 8m telescopes, the HET 
9.2m telescope, the SALT 9.2m telescope and the twin 
Keck 10m telescopes. Efforts are underway to build even 
larger 20m-30m class telescopes: GMT, TMT and ELT. 
Spatial coincidence with a nearby Galaxy will distill the 
large number of counterpart candidates to a small num- 
ber that c an be promptly followed up spectroscopically 
dKulkarni fc Kasli wal 2009). 

lb illustrate the diversity of follow-up scenarios, we 
consider below five case studies of NS-NS mergers. In 
each case, we discuss optimal strategies for identifying 



the EM counterpart of the NS binary merger. Finally, 
we discuss how we can leverage volume information to 
aid EM follow-up strategies for a population of NS binary 
mergers. 

6.1. Individual binaries 

We first examine sky localization and volume errors 
for one beamed NS-NS binary merger at 391 Mpc, and 
four NS-NS mergers that have distances less than 200 
Mpc and lie within the CLU catalog used in this work. 
We choose the five NS-NS mergers described below be- 
cause their geometric properties or sky locations repre- 
sent useful bounds that illustrate the challenges for any 
EM follow-up. The five case studies comprise NS-NS 
mergers with: i) its orbital angular momentum vector 
face-on towards the Earth, ii) a close-by event, iii) a 
source position at low Galactic latitude, iv) a source po- 
sition at high Galactic latitude, and v) a source position 
in a dense galaxy cluster environment. 

6.1.1. Case Study I: Beamed binary merger at 391 Mpc 

We consider the case of a binary merger beamed to- 
wards us. Given the Malmquist bias (§4.3), these bina- 
ries are at threshold and are thus, on average, located 
further away. Out of 200 randomly sampled mergers de- 
tected with a GW Net3a, the distances of beamed NS-NS 
mergers are 391 Mpc, 506 Mpc, 560 Mpc and 564 Mpc. 
Illustrated by Figure [9) using GW networks 3 and 5, the 
localization for the closest of these binaries is 483 deg^ 
(95% c.r.) and 13 deg^ (95% c.r.) respectively. Cur- 
rent 7- and X-ray satellites are easily sensitive to SGRBs 
at these distances (the furthest detected Swift SGRB is 
090426 at a redshift of 2.68 or 22 Gpc). Advantageously, 
these satellites have large instantaneous FOVs. More- 
over, given the precise timing of the ^a mma ray burst, 
false- positive signals are not a concern (Kanner et al. 
2012[ ). If a precise position (e.g. with XRl' onboard the 
Swift satellite or with MIRAX-HXI) is available, prompt 
follow-up to look for the radio and optical afterglows 
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Fig. 7. — The top panel illustrates the depth versus area trade- 
off for optical telescopes when searching for NS-NS mergers de- 
tected by a GW NetSb, where the isotropic optical counterpart 
is assumed to be brighter than Mj^= —14 for at least two hours. 
Bottom panel illustrates the depth versus area trade-off for op- 
tical telescopes when searching for NS-NS mergers detected by a 
GW NetSb and where the optical EM counterpart is assumed to 
be brighter than Mr = — 11 for at least two hours. The different 
colors represent different telescope apertures: Red is l-3m class. 
Green is 4-7m class. Blue is 8-lOm class telescopes. Open square is 
a small FoV camera and filled circle is a large FoV camera in that 
aperture class. Specifically: red square is PTF, red circle is ZTF, 
green square is CFHT, green circle is DECAM, blue square is HSC 
and Blue circle is LSST. 

(which win be much brighter than a kilonova signal but 
decay as a power law in time) will be tractable. 

Unequivocally, wide-field 7- and X-ray satellites (e.g., 
Fermi, Swift, Lobster, MAXI, MIRAX-HXI) are cur- 
rently the most promising wavelengths to search for 
EM counterparts of beamed NS-NS and NS-BH mergers. 
However, as Table 1 shows, the beamed NS-NS mergers 
are a very small fraction (~1.5% to ~3%) of the total 
GW-detected population. Hence, coincident GW and 
EM observations of beamed NS binary mergers will be 
rare. 

6.1.2. Case study II: Close in binary at 70Mpc 

Nearby (< 100 Mpc) NS binary mergers provide ex- 
cellent laboratories in which to study strong-field gravity 
astrophysical processes using joint GW and EM obser- 
vations. Such a "golden binary" should result in high 
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Fig. 8. — Relative fraction of detectable isotropic optical counter- 
parts to NS-NS mergers for GW Net3a and NetSb. Note that even 
the small aperture wide field telescopes are sensitive to a significant 
fraction, and DES and HSC are almost as sensitive as LSST. 
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Fig. 9. — Sky localization errors at 9S% cr. using different GW 
networks for a NS-NS binary merger located at ~ 390 Mpc whose 
orbital angular momentum is directly face-on to the Earth. 
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Fig. 10. — Sky localization errors at 95% c.r. for a NS-NS merger 
located at 69 Mpc with an inclination angle = 150° observed by 
different GW networks. The expected SNRs at LIGO Hanford, 
LIGO Livingston, Virgo, LIGO India and KAGRA are 35, 47, 26, 
24, 36 respectively. 



SNR detections in GWs and multiwavelength EM waves, 
enabling excellent characterization of the physical prop- 
erties of the progenitor and post-merger remnants. 

Let us consider a simulated NS-NS merger located at 
69 Mpc. A GW network 3 measures the source posi- 
tion to 0.6 deg^. A GW network 5 reduces this sky area 
error by a factor of a couple to 0.3 deg^ (Figure 10). Us- 
ing GW network 1, the distance range is from ^Tto 73 
Mpc at 95% c.r. With relatively small localization errors 
and distance measures, the number of astrophysical false- 
positive events that require classification will be nearly 
zero. In addition, assuming that the NS binary merger 
occurred near or within a galaxy, cross-correlating GW 
localization errors with galaxy catalogs, such as the Cen- 
sus of the Local Universe (CLU), leaves us with only a 
han dful of candidates for galaxy hosts. As Figures 11(a) 



and 



11(b) 



show, 



we expect to see five galaxies within 
an error cube of 2 degx2.5 degx55 Mpc. The CLU is 
currently 65% complete within this distance bin and ef- 
forts are underway to ma ke this catalog more complete. 
With GW network 5, Fig. |1 1(b)] shows that we should be 
able to identify uniquely one host galaxy using full 3D- 
marginalized volumes computed at 95% c.r. (this is also 
the case using the low-latency GW volumes described in 

Therefore, for such a well-localized nearby "golden" 
binary (with distances < 100 Mpc), we can undertake 
pointed host galaxy follow-up and are no longer limited 
to large FoV cameras. Moreover, we can observe deeper 
and use a faster cadence than for an average faint binary. 
Intensive panchromatic follow-up with a wide array of 
facilities in the optical, infrared, radio and X-ray would 
ensure that we leave no stone unturned in studying the 
EM counterpart to a golden binary. We note that only 
10 % of NS-NS mergers seen by GW Net3a are golden 
and will have true distances less than 100 Mpc. Using 
GW Net5b, we note that only 1.5% of binaries will be 
golden, but they will have high SNR detections in GWs 
and EM waves. 
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Fig. 11. — Volume errors at 95% c.r. for a NS-NS merger located 
at 69 Mpc with inclination angle = 150 ° observed by GW networks 
3 and5. The blue circle marks the merger's true position and the 
black circles denote the 5 galaxies within the cube's volume. Note 
that the above circles are not to scale. 

6.1.3. Case Study III: Low Galactic Latitude and at 125 

Mpc 

Let us consider the case of a binary with inclination 
angle 150°, distance of 125 Mpc and located very low on 
the Galactic plane (latitude of -0.11°, longitude of 63.4°). 
Using GW network 3, the sky localization of the merger 
is 1.8 deg^ at 95% c.r. GW network 5 reduces the sky 
area error to 1.3 deg^ . 

Without any distance information, we would need to 
search for the isotropic Mi^=— 14 optical counterpart sig- 
nal out to 450 Mpc, i.e. 24.3 apparent mag. The number 
of background supernovae active in this area at this time 
would be ~ 120 and the number of foreground M-dwarf 
flares active would be 76. Unfortunately, at such a low 
Galactic latitude, galaxy catalogs are most incomplete 
and cannot be effectively used to reduce false positives. 

However, the derived GW low-latency localization vol- 
ume can help reduce false positives. Using GW network 
3, the distance measure in this volume ranges from 72 
Mpc to 142 Mpc at 95% c.r. (Figure [T2|. The reduc- 
tion in volume from the maximum detectable distance of 
450 Mpc to 142 Mpc is 97%! Instead of searching down 
to 24.3 mag, we only need to search to 21.8 mag. This 
reduces the background and foreground false positives to 
< 10. 

Another effective strategy to deal with foreground false 
positives is to use a quiescent star catalog (in the optical 
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Fig. 12. — Volume errors at 95% c.r. for a NS-NS merger lo- 
cated at 125 Mpc with inclination angle = 150° observed by GW 
networks 3 and 5. The blue circle (not to scale) marks the bi- 
nary's true position. Due to the low Galactic latitude, there are 
no known galaxies within the volume shown. Using network 3, the 
GW distance measure ranges from 72 Mpc to 142 Mpc (95% c.r.); 
the reduction in volume is 97%. 

or infrared) that is about two magnitudes deeper than 
the search depth (Stubbs 2008). Several synoptic surveys 
(e.g., SDSS, PTF, PanStarrs, Skymapper in the optical; 
VISTA, UKIDSS, WISE in the infrared) are underway to 
give us such a catalog. Moreover, some of these surveys 
will also provide a multi-year historic baseline for variable 
sources. 

In practice, the limiting factor for follow-up of such 
a merger would be the crowding and large line of sight 
extinction; ~ 8 apparent mag in i?-band and 1 apparent 
mag in K-band. Near-infrared follow-up would be much 
easier than optical follow-up. The percentage of GW 
detected mergers with a Galactic latitude less than 5° is 
~ 9% and Galactic latitude less than 10° is ~ 18%. 

6.1.4. Case Study IV: High Galactic Latitude and at 139 

Mpc 

The previous three case studies represent a subset of 
binaries. Let us now consider a canonical binary, which 
is neither beamed, nor very close-by, nor too low on the 
Galactic plane. This binary has an inclination angle of 
64°, a distance of 139 Mpc and a Galactic latitude of 
—66°. We expect no EM counterparts at 7- and X-ray 
wavelengths. Using GW network 3, GW measurements 
can localize the event to 19.5 deg^ on the sky (using GW 
network 5, the localization error improves to 8 deg^). 

Without any GW distance constraints, we would need 
to search for an optical counterpart brighter than Mr 



= —14 to a horizon distance of 450 Mpc i.e. 24.3 ap- 
parent mag. Therefore, in a 19.5 deg^ error circle, the 
extragalactic false-positive number will be ~ 1300. The 
Galactic false-positive number will be 100. 

One strategy is to identify the EM counterpart from 
a false positive based on the light curve signature. An 
outburst due to a NS binary merger would be a one-time 
occurrence. If we have a good historic light curve of the 
candidates, foreground false positives and AGN would 
show previous eruptions. Unfortunately, at a depth of 24 
apparent mag, this may not be the case for most of the 
sky until LSST is operating for a few years. We can also 
use the theoretical prediction that optical counterparts, 
such as kilonovae, evolve faster than supernovae in the 
same field. Unfortunately, if we wait too long to obtain 
multiple epochs, the EM counterpart may fade to a level 
where it is too faint for spectroscopic follow-up. We could 
also use theoretical predictions that kilonovae may be 
redder than supernovae in the same field. Unfortunately, 
given the depth needed and large localization areas, there 
may not be enough time to obtain data in multiple filters. 

To reduce the number of false positives, a simple ap- 
proach is to assume that the merger is spatially coinci- 
dent with or nearby a galaxy within the distance reach 
of the GW network. Galaxies occupy a very small area 
on the sky. There are only 228 known galaxies in the 
error circle within 450 Mpc. Allowing a large radius in- 
cluding 50kpc around each galaxy, the total area is 60 
arcmin^ and the reduction in false positives is a factor of 
1200! Allowing an even larger radius including 100 kpc 
around each galaxy, the total area is 240 arcmin^ and the 
reduction in false positives is still a factor of 300. Unfor- 
tunately, the current galaxy catalog is grossly incomplete 
and it is imperative that efforts be made to complete it. 

Furthermore, we can leverage our localization volume 
constraint to further reduce the number of relevant host 
galaxies. Using GW network 3, the merger's distance 
is measured to be between 108 Mpc and 228 Mpc, the 
volume is smaller than between OMpc and 450 Mpc by 
88%. A smaller volume would correspond to a smaller 
number of galaxies. In this case, 73 galaxies are known 
to lie in the localization volume. Using GW network 5, 
the merger's distance is measured to be between 117 Mpc 
and 230 Mpc and currently, ~ 37 galaxies are known in 
this volume. This list of galaxies is incomplete by a factor 
of two. If we had a complete catalog of galaxies, and for 
cases where the number of galaxies was less than few tens, 
we would even consider targeting galaxies individually. 
We would not require a large FoV camera, paving the 
way for using infrared, radio and X-ray facilities. 

6.1.5. Case Study V: Dense galaxy cluster environment at 
high Galactic latitude at 115 Mpc 

For our final case study, we choose a simulated NS-NS 
binary merger event that occurs in an extremely dense 
galaxy cluster (ABELL 4038). Here, the use of localiza- 
tion volumes to target individual galaxies is not feasible 
as the number of galaxies is too large. However, using 
a galaxy catalog to prioritize follow-up can still reduce 
false positives by orders of magnitude. 

This specific binary has a distance of 115 Mpc and 
Galactic latitude and longitude of —75.8° and 25.7° re- 
spectively. With GW network 3, its sky localization can 
be measured to 18.8 deg^ and its distance measure ranges 
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Fig. 13. — Volume errors at 95% confidence interval for a NS 
binary merger located at 139 Mpc with inclination angle = 64° 
observed by GW networks 3 and 5. The blue circle marks the 
binary's true position and the black circles denote the 279 galaxies 
within the cube's volume. Note that the circles above are not to 
scale. 

from 110 Mpc to 269 Mpc. The upper measure of the dis- 
tance (at 95% c.r.) is greater than the maximum distance 
of 200 Mpc used in CLU. Assuming only a horizon dis- 
tance of 200 Mpc and using the CLU catalog, we find 780 
galaxies comprising an area of 1070 arcmin^ within such 
an area of the sky. If we then include our lower measure 
of the distance (at 95% c.r.), we find 410 galaxies com- 
prising an area of 560 arcmin^. Increasing the network 
from three to five interferometers, its sky localization im- 
proves to 14.4 deg^ and its distance measure ranges from 
97 Mpc to 155 Mpc. Using CLU, we then find 390 galax- 
ies comprising an area of 530 arcmin^ within such an area 
of the sky. 

In such a dense galactic environment, targeting hun- 
dreds of galaxies individually is not feasible. However, 
the reduction of false positives is still significant. Specif- 
ically, a snapshot of 18.8 deg^ out to 450 Mpc, would give 
1250 background supernovae. Restricting the search to 
candidates spatially coincident with nearby galaxies (i.e. 
1070 arcmin^) reduces the false positives to 20. Further 
imposing the volume constraint of 269 Mpc, reduces the 
false positives to 4 events. 

6.2. Astrophysical populations 

As discussed in §6.1[ the combination of GW distance 
information, GW sky localization and galaxy catalogs 
could help reduce the number of galaxies that are possi- 
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Fig. 14. — Volume errors at 95% confidence interval for a NS- 
NS binary merger located at 115 Mpc with inclination angle = 
87° observed by GW networks 3 and 5. The blue circle marks the 
binary's true position and the black circles denote the 1350 galaxies 
within the cube's volume. Note that the circles above are not to 
scale. 



ble hosts of the NS binary merger event and the number 
of false-positives. The use of spatial coincidence with 
galaxies in turn could have a substantial effect in reduc- 
ing the number of false-positive transients that ne ed to be 
considered as possibl e NS-binary merger events (Kulka- 
rni fc Kasliwal||2009 ) . Alternatively, if we can limit the 
number of host galaxy candidates to a few as in the case 
of a golden nearby GW merger, targeted follow-up op- 
portunities for individual galaxies becomes a possibility. 

next, we com- 



Following from Figures 6(a) and 6(b) 



pare how well we fare when computing GW volumes ei- 
ther using the event's distance measure or using only the 
GW horizon distance. In particular, we compute TvoI the 
fractional change in volume: 



(min[dn,d/i])^ 



d^ 



(8) 



for samples of NS binary mergers detected using different 
GW networks and triggering criteria (see Figures 15(a) 



and [15(b)] ) 

We find that the fractional change correlates with dis- 
tance to the binary and the reduction is higher for nearer 
binaries. Specifically, for NS-NS mergers, the fractional 
change in volume is < 60% for those events located 
within 200 Mpc (-- 1/3 of the detected binaries). For NS- 
BH binaries, the fractional change in volume is < 60% 
for those events located within 700 Mpc (~ 1/2 of the 
detected binaries). 
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Fig. 15. — Fractional reductions in volume for a detected sample 
of NS-NS mergers observed by GW Net3a (top panel) and NS-BH 
mergers observed by GW Net 5b (bottom panel). This fraction is 
the ratio of the volume encompassed by GW distance measures 
(Figure 6 in Section 4.5) and the volume encompassed within the 
GW maximum detectable distance. Each circle represents a de- 
tected NS binary merger. The different colo rs re present different 
case studies of NS-NS mergers examined in §6.11 green square is 
Case I (beamed binary), red circle is Case 11 (a nearby binary), 
light blue triangle is Case 111 (a merger at low Galactic latitude), 
purple diamond is Case IV (a merger at high Galactic latitude), 
and orange star is Case V (a binary in a dense galaxy cluster en- 
vironment). 

The implications for search strategies are as fohows. 
In cases where the upper hmit of locahzation volume is 
much less than GW maximum detectable distance (c.f. 
case studies 2, 3 and 4), the EM counterpart search can 
be to brighter apparent magnitudes. This significantly 
reduces false positive numbers and spectroscopic follow- 
up is easier. More frequent are the cases where knowledge 
of the lower distance limit reduces the total number of 
host galaxies. 

The number density of galaxies is 0.6 per deg^ per Mpc 
at z= Q.l and L>0.1 and scales with (Blanton et al. 
|2QQ3p . For a median localization of 10 deg^ and out to 
a distance of 200 Mpc, the number of galaxies is ~80. 



The median size of a galaxy is [0.8, 0.4, 0.2] arcmin^ 
at [50, 100, 200] Mpc respectively. The diameter as- 
sumed here is for the surface brightness contour of 25 
mag arcsec"^. This corresponds to a projected offset 
(offsetkpc=sizerad X distancckpc) of [7.3, 10.8, 14.0] kpc 
at [50, 100, 200] Mpc respectively. Thus, the total area 
occupied by 80 galaxies is 0.004 deg^, a factor of 2500 
smaller than the localization. We can easily search 10 
times the size of the galaxy to accommodate largekicks 
of the order of h undred kiloparsec (e.g., Fong et al.|2010 



Kasliwal et aT1|2012J in NS-NS mergers and still gain a 



factor of 25 in terms of reduction of false pos itives. How- 
if the kicks are over a megaparsec (Kelley et al. 



ever, 



2010), we cannot use the positions of host galaxies to 
reduce false positives. 

Efforts are underway to complete the CLU galaxy cata- 
log using four narrow-band filters on the Palomar 48-inch 
Schmidt telescope. This will boost the completeness from 
50% to 85% of the 5-band light at 200 Mpc in the three- 
quarters of the sky accessible from Mount Palomar. Note 
that this distance limit is well-matched for majority of 
NS-NS binaries detected by GW Net 3a. However, par- 
ticularly for NS-BH binaries, we should consider an even 
larger effort to complete galaxy catalogs out to several 
hundred Mpc. 

7. CONCLUSIONS AND DISCUSSION 

Observing compact-object binary mergers in both GW 
and EM will be challenging. GW interferometers will 
only be able to localize the merger to sky errors ranging 
from tens to hundreds deg^. In addition, the estimated 
rates of GW-detected compact binary mergers span sev- 
eral orders of magnitude from zero to hundreds of merg- 
ers per year. Theoretical predictions of EM signatures 
from the optical to radio is also an active area of re- 
search and model estimates vary significantly. Finally, 
our characterization and understanding of the transient 
sky at different wavelengths, timescales and sensitivity 
remains incomplete. 

In this paper, we present the first comprehensive 
end-to-end simulation of the detectability and identi- 
fication of neutron star binary mergers with GW and 
EM facilities. Our simulation comprises: the construc- 
tion of astrophysically-distributed populations of GW- 
detectable NS-NS and NS-5 M© BH binary mergers, GW 
source characterization using different GW detector net- 
works and triggering criteria, establishing the detectabil- 
ity of plausible EM counterparts by upcoming or cur- 
rent telescopes (particularly, optical), and identifying the 
GW event among the few-to-many astrophysical false- 
positive transients in different wavelengths. The extent 
of our analysis is naturally dependent on the assumptions 
made: an underlying population of non-spinning NS-NS 
or NS-5 M0 BH binary mergers, GW instrumental Gaus- 
sian noise with each interferometer operating continu- 
ously, negligible GW systematic errors arising from the 
GW waveform and instrument calibration, a high avail- 
ability factor for telescopes to follow-up GW mergers, 
and idealized optical observing conditions. 

Our work is novel in six principal ways. First, we 
construct GW-detected populations of astrophysically- 
distributed NS binary mergers using CLU (a volume- 
limited local universe galaxy catalog out to 200 Mpc), 
instead of at fixed SNR or distances as in earlier works. 
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Second, we incorporate advanced LIGO detectors us- 
ing optical squeezed light into our analysis. Third, we 
compute explicit marginalized 3D GW volumes using 
MCMC. Fourth, we consider how GW volumes can assist 
follow-up by optical telescopes and other EM facilities. 
Fifth, we quantify the tradeoff between depth and area 
for a variety of optical telescopes, including 4- 7m class 
and HSC telescopes for the first time. Sixth, by exam- 
ining individual NS binary mergers, we suggest how to 
pinpoint the GW event amongst a possible plethora of 
astrophysical Galactic and extragalactic false-positives. 

By expanding the parameter space of NS binary popu- 
lations and GW networks, we show that GW detectable 
distances and sky area errors may range an order, or 
several orders of magnitudes respectively. From our case 
studies of GW-detected NS binary mergers, we find that: 

• Thanks to the GW Malmquist effect, the fraction 
of NS binary mergers beamed towards us (with 
Oj < 6°) is boosted to large distances (400 Mpc-1.3 
Gpc). However, we have shown empirically that 
the fraction is still tiny (1-2.5%). For this subset 
of events, the easiest identifiable EM counterpart 
would be a contemporaneous SGRB. All-sky 7-ray 
detectors are essential to ensure joint GW and EM 
observations. Optical squeezing implemented in 
GW interferometers will also increase the beamed 
fraction by a further factor of nine to ten. 

• For the small number of golden nearby binaries 
(< 100 Mpc), which have small GW localizations 
given the high SNRs of GWs, we should inten- 
sively follow-up such events at all wavelengths. 
Whilst the number of detected golden binaries is 
independent of network, the total GW-detected 
fraction ranges from 2%-10% from GW NetSa to 
GW NetSb. The number of false-positives in their 
small localization areas or volumes should be small 
(§6.1.2). 

• For the majority of binaries which are neither 
beamed nor nearby, the challenges in detecting the 
isotropic optical counterpart are surmountable by 
optimizing the depth versus area tradeoff. Initially, 
we expect small telescopes (especially given their 
larger number and wide-field cameras) to play an 
important role in detecting counterparts. GW net- 
works with fewer detectors will have poorer local- 
izations and lower distance sensitivity. Small tele- 
scopes have strength in numbers and can be ex- 
pected to be more flexible for rapid follow-up for 
GW triggers. Later, with increasing detector num- 
bers and instrument sensitivity and use of coherent 
triggers, GW networks will have higher distance 
sensitivity and improved localizations, larger tele- 
scopes will be essential (§5.2.1). Given weather and 
limited sky accessibility of optical telescopes (due 
to sun constraints and altitude constraints) , we ad- 
vocate for a world-wide network of telescopes of dif- 
ferent sizes with wide-field cameras spread across 
different latitudes and longitudes. Furthermore, we 
advocate the building of larger FoV infrared cam- 
eras. 

Our simulation and detailed case study analysis moti- 



vate us to search for EM counterparts to GW binaries. 
We find that although there are challenges, they are sur- 
mountable by timely advance preparation. Hence, we 
conclude here with four action items that will better pre- 
pare us to securely identify the detected EM counterpart: 

1. We should complete as much as possible host 
galaxy catalogs out to z ~ 0.1 to increase iden- 
tifying an EM counterpart through two different 
means: (i) spatial coincidence with a nearby galaxy 
can quickly eliminate false positives for a subset of 
NS binary mergers by orders of magnitude. This is 
critical to prioritize candidates for prompt spectro- 
scopic and panchromatic follow-up (§6.1.4, §6.1.5, 
§6.2). (ii) together with low-latency GW volume 
errors, in some cases, the number of galaxies can 
be reduced to a tractable number for possible tar- 
geted follow-up with relatively smaller field-of-view 
facilities e.g. radio, infrared. X-ray, large aperture 
optical telescopes (§6.1.2, §6.1.4, §6.2). An ongo- 
ing effort is the Census of the Local Universe (CLU) 
using narrow-band filters (Ha) on the Palomar 48- 
inch. Another planned survey (H I) is WA LLABY 
with ASKAP in the Southern Hemisphere (jKaplan 



2012). 



2. We should construct deep (~ 26 apparent mag) all- 
sky quiescent stellar source catalogs which would 
help eliminate foreground false positives, particu- 
larly at low Galactic latitudes (§6.1.3). We esti- 
mate that we will require catalogs that are approx- 
imately two magnitudes deeper than the EM coun- 
terpart and span optical and infrared wavelengths. 
Current ongoing efforts may not be deep enough, 
e.g. SDSS, PTE, PanStarrs, SkyMapper, WISE, 
VISTA. 

3. We should rehearse the search for transients in 
large sky localizations. For example, the successful 
identification and follow-up of an optical afterglow 
of a Fermi/GBM G RB in a three de^.^ error circle 
(c.f. GRB120716A, |Cenko et "ar] r2012) is encour- 
aging. Such efforts are a full dress rehearsal for 
elusive EM counte rparts in GW co nstrained large 
swaths of the sky ( Kulkarn"i||2012a ) . 



4. We should continue to construct a complete inven- 
tory of transients within several hundred Mpc. Just 
in the past few years, we have uncovered multiple, 
new classes of optical transients, which are fainter, 
faster and rarer than supernovae (e.g., jKasliwal] 
2011|) . We may indeed even be lucky enough to see 



an EM counterpart to a NS binary merger prior to 
hearing the GWs! 

In summary, given the diversity of properties and lo- 
cations of possible EM counterparts and challenges in 
their identification, we advocate a comprehensive, multi- 
wavelength multi-pronged approach to observe compact 
binary mergers in GWs and EM waves. 
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